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ABSTRACT 


An analytical and experimental investigation of two 
phase flow in a discharge line with a sharp edged entrance 
Xj presented for both vibratory and non-vibratory conditions. 

Initial efforts were devoted to single phase liquid flow 
in the discharge line. Classical one-dimensional fluid flow 
relations yielded analytical predictions which were closely 
matched by experimental data. While no attempt was made to 
develop predictions in the vibratory cases, an extensive 
experimental investigation showed that vibration of the test 
section generally caused an increase in pressure losses. 

A modified cavitation number was developed based on a 
free streamline analysis which considered the flow into the 
sharp-edged entrance as flow from an orifice. Good agreement 
with experimental data was sho wn to exist for several differ- 
ent geometric configurations. Furthermore, a time dependent 
cavitation number was developed for the cases where vibration 
was present. An analytical model describing choked two-phase 
flow was developed using the free streamline theory. Excel- 
lent agreement between theory and experimental data was 
achieved. No effect of vibration on choked flow was detected 
experimentally. Expressions associated with the normal shock 
separating the two-phase fluid (liquid-vapor) and single 
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phase fluid (liquid) were derived. In particular, equations 
defining the minimum pressure drop for choked flow and equa- 
tions predicting the shock position as a function of the 
stagnation pressure loss are given. Good agreement between 
the analytical predictions and experimental data was obtained 
in both cases. 
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CHAPTER 1 
INTRODUCTION 

High speed transfer of liquids is required for many 
modern engineering applications. Numerous examples can be 
cited in rocket propulsion, automatic food processing manu- 
facturing, rapid emergency shut-down of nuclear reactor power 
plants, and in several other areas of technology. 'In the 
development of equipment required for these applications, 
three major fluid flow problems have been encountered. First, 
rapid liquid transfer has often been accompanied by unaccept- 
ably large frictional losses, as would be expected. Second, 
local regions of very high velocity have induced liquid cavi- 
tation resulting in adverse operating characteristics, and 
even causing ultimate damage to the system in some instances. 
Third, totally two-phase flows have appeared in situations 
where they were not expected, usually producing very detri- 
mental effects which severely limit system performance. The 
last two problems have been particularly important for cryo- 
genic liquids where heat transfer typically contributes to 
significantly increasing the amounts of vapor formation. In 
some cases vibration of equipment due to external sources 
(sxich as combustion instabilities in rocket motors or im- 
properly balanced rotating machinery) has affected fluid flow 



2 


rates significantly, but in ways which previously have not 
been well understood. 

The research described in this thesis was concerned 
with fluid flow in tubes with a view toward a study of flow 
conditions occurring in propellant tank discharge lines of 
liquid fuel rockets. The overall investigation comprised 
three major topics: 

1. Single phase liquid flow 

2. Cavitation inception (onset of two-phase flow) 

3. Totally two-phase flow accompanied by liquid 
flashing, choking, and normal shock phenomena 

Each of these flow regimes was studied both in the absence 
and in the presence of longitudinal vibration of the system. 

For accomplishing the goals of the research in each of 
the three areas listed above, an experimental facility was 
designed and constructed. Major initial efforts were devoted 
to a study of its operating characteristics and the function- 
ing of its associated instrumentation. The experimental 
portion of Topic 1 (single phase liquid flow) was then begun. 
The classical predictions for single phase liquid pipe flow, 
flow through an orifice, and floi* through a restriction were 
used as the established information against which the experi- 
mental data were compared. The details of this initial single 
phase flow research are described in Chapter 6. 

The second major portion of the investigation was con- 
cerned with the generation of two-phase flow as the result 
of cavitation. The classical cavitation parameter was used 
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as a basis from which a more comprehensive relation was de- 
veloped for correlating cavitation inception. A modified 
cavitation number was also formulated to account for cavita- 
tion in the presence of time dependent pressure oscillations 
resulting from vibration of the system. The analytical 
results of the cavitation study are contained in Chapter 7 
along with the appropriate experimental verifications. 

The third and final major portion of the investigation 
was associated with the experimental observation of flashing 
flow phenomena and the corresponding analytical predictions. 
Two unique correlations were developed which were used to 
predict the occurrence of choked flow as well as the incep- 
tion and position of a normal shock. These items are dis- 
cussed in Chapter 8. 

To summarize the various topics outlined above, a dia- 
gram illustrating the organizational arrangement of the 
entire investigation is given in Figure 1. While the results 
of this research (Chapters 6-8) are pertinent to many flow 
applications, the initial interest was prompted by the 
problems associated with liquid fuel rocket systems. .For this 
reason, Chapter 2 contains a general description of several 
two-phase flow phenomena, and Chapter 3 contains a description 
of the background of some typical fluid flow problems associ- 
ated with liquid propellant rocket technology. These chapters 
are followed by a literature survey in Chapter 4 and a descrip 
tion of the experimental facility in Chapter 5 prior to the 
presentation of the three major topics of the present research 
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CHAPTER 2 

GENERAL DESCRIPTION OF LIQUID FLOW 
AND LIQUID-GAS FLOW PHENOMENA 

In this chapter liquid flow and liquid-gas flow are 
described, and some of the useful analytical methods for 
treating them are briefly discussed. Consider a liquid 
flowing through a constriction in a line or past a sharp 
edged boundary. A localized region of high velocity and 
low pressure is present in either situation because of the 
necessary forces acting on the liquid to produce its rapid 
acceleration in such a region. If the pressure becomes low 
enough, the liquid may vaporize causing the onset of two- 
phase flow. Once two-phase flow occurs in these low pressure 
regions, the resulting bulk flow rate of the fluid can be 
changed very drastically in some instances. 

A similarity analysis of adiabatic cavitation, which 
involves vapor generation in the absence of heat transfer, 
has led to the formulation of the cavitation number, a. This 
parameter can be used to predict the onset of two-phase con- 
ditions in a given system; however, its effectiveness in such 
an application is limited by the condition of the liquid. 
Dissolved gases and entrained gas bubbles, as well as solid 
particles suspended within the liquid, tend to affect the 
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inception of cavitation. Consequently, these factors lessen 
the usefulness of the cavitation number in accurately pre- 
dicting the onset of two-phase flow. 

In high speed liquid flows, cavitation and other forms 
of two-phase generation can severely limit the validity of 
a single phase liquid flow analysis. The major difficulty 
lies in determining the limits of application of a predeter- 
mined liquid flow solution. Whether analytical solutions 
for the flow patterns are obtained from the use of Bernoulli's 
equation for incompressible flow, conformal transformation 
techniques, or very lengthy numerical solutions to the Navier- 
Stokes equations, they necessarily become invalid at the 
onset of a two-phase flow. 

Since a high speed liquid flow can generate a two-phase 
flow situation, it also becomes desirable to obtain analyti- 
cal descriptions of these more complex two-phase flow pheno.-.- 
mena. However, the analysis of two-phase flow has been found 
to be quite difficult in general. Only recently has a con- 
centrated effort been directed towards obtaining solutions 
to basic relations for two-phase flow. Many models have been 
postulated to describe the different flow regimes that characr 
terize two-phase systems and to compare them with equivalent 
single phase systems. A few of the generally accepted flow 
patterns are classified as homogeneous, bubbly, stratified, 
annular, wavy, slug, plug, and spray. These patterns exist 
separately or in various combinations (such as spray-annular 
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or .-homogeneous bubbly flow) , and the patterns may exist 
in different forms for vertical and horizontal flow situa- 
tions. Generally speaking, the interactions between the two 
phases are not completely understood. Therefore, the rela- 
tions for two-phase flow are not considered to be as well 
substantiated as those for single phase flow. 

Choking of two-phase flows can occur just as in the 
case of single phase flow, but the resulting maximum two- 
phase flow rate is in general much lower than the maximum 
flow rate for the corresponding single phase flow. This 
phenomenon occurs when the fluid velocity equals its acoustic 
velocity at some location in the flow path. For an ideal 
incompressible liquid, the acoustic velocity is theoretically 
infinite. However, any real liquid that is typically assumed 
incompressible in reality does have some compressibility, 
and therefore possesses a correspondingly large, but finite 
value of the acoustic velocity. When a gas or vapor is 
present in a nearly incompressible liquid, the resulting 
two-phase fluid mixture displays some of the physical be- 
havior of a highly compressible fluid. The available analy- 
tical expressions which predict choking conditions for such 
a two-phase fluid yield different results depending upon the 
type of flow pattern that is assumed. The homogeneous flow 
model is the most direct to use in deriving an acoustic 
velocity relation. This model is based on the concept of a 
pseudo fluid having homogeneous properties and a single, 
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definable txtfo-phase acoustic velocity. In this formulation 
the two phases are combined, and the mixture is treated as 
a single, compressible fluid. The other models usually con- 
sider the acoustic velocities in the two phases as being 
different from each other. For the homogeneous mixture 
model, the predicted acoustic velocity exhibits a functional 
dependency on the void fraction and can be shown to have a 
minimum value at a void fraction of one-half. 

When a saturated liquid flov/s into a region where the 
pressure is lower than the vapor pressure of the liquid, 
some of the liquid can "flash” into vapor resulting in a 
two-phase mixture. This phenomenon occurs due to very rapid 
property changes of the fluid which are brought about by 
bulk nucleation within the liquid. The validity of the 
homogeneous model for predicting the acoustic velocity and 
the resulting choking conditions becomes questionable when 
the homogeneous mixture results from flashing. The assump- 
tion that the vapor phase obeys the relations which describe 
compressible fluids is questionable because of the non-equi- 
librium effects of the flashing process. 

The lack of a proper definition for the acoustic veloc- 
ity in a flashing flow severely limits the usefulness of 
the acoustic velocity criterion for predicting the corres- 
ponding choking conditions analytically. A general defini- 
tion for choking to be used in the remainder of this thesis. 



and which can be applied to any type of flow, is the fol- 
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A flow system is choked when the mass flow 
rate reaches a maximum value for given upstream 
stagnation conditions which is not exceeded re- 
gardless of the downstream conditions. This 
implies that there is a critical exit pressure 
at which the maximum flow rate is attained, and 
a decrease in exit pressure below this critical 
value causes no further increase in the flow 
rate for the given, fixed, upstream stagnation 
conditions . 

Obviously this definition can be applied to a gas, liquid, 
or a two-phase fluid, including one that is flashing, even 
though the non-equilibrium conditions in a flashing flow are 
not fully understood. 

While the flow phenomena such as cavitation, flashing, 
and choking can be described on an individual basis to a 
limited degree of accuracy, the application of this know- 
ledge should be extended to situations where interdependen- 
cies between them exist within a single flow system. The 
threshold conditions accompanying the onset of a two-phase 
mixture in an otherwise pure liquid flow should be investi- 
gated and specified for use in those areas of technology 
where these situations develop. That is, the limiting con- 
ditions of single phase liquid flow should be determined, 
and the initial conditions within the two-phase mixture 
should be specified to provide the necessary correlation 
between the two types of flow. Also, for any given flow 
system, the relationship between the two-phase phenomena 
and their overall effect on the flow and the overall system 



10 


performance should be described. That is, the pressure 
losses, flow rates, and velocities should be related to the 
many different flow conditions that can exist within a flow 
line. Only in this manner can effective engineering design 
of these systems be accomplished. 

The following chapter describes the relationships be- 
tween an application area of considerable interest (liquid 
fuel rocket technology), high speed liquid flow phenomena , 
and the research described in this investigation. 
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CHAPTER 3 

RELATION OF RESEARCH 
TO LIQUID FUEL ROCKET TECHNOLOGY 

The effects associated with high velocity liquid flows 
are ever present in liquid fuel rocket technology. This is 
due to the low temperatures of the cryogenic liquids used 
and the resulting high heat transfer rates that typically 
occur. In the case of the onset of two-phase flow, vapori- 
zation is brought about when the liquid reaches its saturation 
temperature corresponding to the local pressure. Thus, the 

presence of heat transfer enhances the occurrence of this 

l 

condition. As a result of heat transfer, vaporization of 
liquid can occur in propellant tanks by boiling, and cavita- 
tion can be induced in the lines, both having possible detri- 
mental effects on the total propellant discharge system per- 
formance . 

Large flow rates are a necessity in some of the advanced 
rocket designs. For example, the first stage of the Saturn V 
rocket consumes approximately fifteen tons of propellants 
during each second of operation. The propulsion systems of 
liquid fuel rockets, including fuel tanks, discharge lines, 
pumps, valves, and associated hardware items, are designed 
to handle the propellants in their liquid phase. When 
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two -phase flow occurs due to vaporization of the liquid, 
serious effects on the performance of the total system can 
result . 

Two-phase conditions which develop in propellant supply 
systems are believed to contribute to the existence of com- 
bustion instabilities in the rocket motors. These instabi- 
lities can cause severe vibration of the entire structure of 
a rocket. The dynamic effects are sometimes directly respon- 
sible for sizeable pressure fluctuations within the liquid 
propellants. The exact manner in which the structural vibra- 
tional characteristics of discharge lines affect liquid 
flow, as well as the inception of two-phase flow, is not well 
known or understood at present. 

Certain dynamic effects of vibration on tanks, partially 
filled with liquid and pressurized with an ullage gas, have 
been investigated and explained previously, at least in a 
qualitative manner. The pressure oscillations within the 
fluid are dependent upon frequency, acceleration level, 
location in the tank or supply line, and structural rigidity. 
The stability of the interface between the liquid and the 
ullage gas in the tank is dependent upon the same parameters. 
Under certain specific conditions, considerable gas entrain- 
ment occurs at this interface resulting in large conglomera- 
tions of gas bubbles which subsequently become located at 
low pressure regions within the fluid. As the liquid level 
in a tank drops during discharge, the dynamic response of 
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the tank changes such that the bubble agglomerations either 
rise to the surface and dissipate into the ullage gas region, 
or are discharged into the supply line to be subsequently 
passed on to a pump, and perhaps also to the combustion 
chamber of the rocket motor. This type of two -phase flow is 
different from the cavitation phenomena described earlier 
since this occurs when the ullage gas is a permanent foreign 
gas which does not undergo a phase transformation to a 
liquid. The resulting liquid-gas flow is more properly 
defined as two -phase, two -component flow. 

The entrainment of ullage gas into the liquid and the 
subsequent tx^o-phase conditions can have a very serious 
effect on the propellent flow rate. The trapped ullage gas 
bubbles act as nuclei for large scale severe cavitation to 
develop. So significant is the effect of these nuclei on 
cavitation that extensive experimental work has been per- 
formed on cavitation resulting from longitudinal vibration 
in liquid-filled containers in the absence of flow. The 
implications of cavitation in a "no-flo\tf" situation are 
significant since cavitation is usually associated with high 
velocity flows. How the effect of longitudinal vibration of 
a simulated discharge line affects cavitation inception 
during flow through the line has not been thoroughly inves- 
tigated previously. 

Experiments have shown that localized flashing can 
occur in a liquid container which is vibrated on an 
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electro -dynamic shaker at various frequencies and accelera- 
tion levels. It has also been shown that gravity driven 
discharge rates from containers exposed to vibration are 
reduced. The subsequent questions raised by these experi- 
ments are: What effect does vibration have on the single 

phase liquid flow through a discharge line? What effect 
does vibration have on choking and two-phase conditions with 
in a discharge line, and how do these conditions affect the 
propellent transfer rates in a rocket system? 

The answers to the above questions necessarily require 
a thorough knowledge of the phenomena which are occurring. 

To determine the current state of knowledge with regard to 
the phenomena involved in liquid fuel rocket technology and 
discussed in this chapter, a literature survey was conducted 
A description of the pertinent contributions in this area is 
given in Chapter 4. 
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CHAPTER 4 
LITERATURE SURVEY 

An extensive research effort has been devoted to the 
field of two-phase flow in the past two decades. Numerous 
papers have been published covering the multitude of problems 
associated with cavitation, acoustic velocity, critical flow, 
boiling, mathematical modeling and many other topics not 
directly associated with this investigation. The references 
in these areas which are pertinent to the present study are 
discussed briefly in this chapter. Phenomena associated with 
a vibrating environment are of a more special nature, and 
these are reviewed subsequent to discussion of literature 
pertaining to two-phase flow. 

Streeter (1)*, Shames (2), and Kenyon (3), describe the 
contributions of many past investigators in the science of 
fluid mechanics, and material from these well known texts 
formed the starting point for the single phase liquid analysis 
in the present study. The classic relations of fluid mechan- 
ics not only were useful for the liquid flow problem, but 
also served as the basis from which extensions to two-phase 
flows were made. 


* Numbers in parenthesis represent references listed in the 

Bibliography. 
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Some of the very early contributions to the two-phase 
flow literature were associated with jet flow analysis of the 
situation in which a non-viscous fluid of large density was 
discharged into a second non-viscous fluid of much lower 
density. A review of work in the jet flow field was present- 
ed by von Mises (4) who also described his own derivation for 
the two dimensional slot floitf solutions he collected. The 
most recent significant contribution to jet flow theory was 
made by Gurevich (5) who has included many references not 
previously available to free world researchers. This text 
also illustrates the use of complex analysis in solving 
certain types of fluid flow problems in two and three dimen- 
sions . 

An experimental investigation by Rouse and Abul-Fetouh 
(6) was very pertinent to the present research.- Their work 
on jet flows for axisymmetric geometries utilized an elec- 
trohydrodynamic analogy. The resulting experimental solution 
to an axisymmetric jet flow problem was shown to be in agree- 
ment with the two dimensional slot flow analysis performed 
by von Mises (4). It was shown that the two different geo- 
metric solutions presented in References (4) and (6] , two- 
dimensional and axisymmetric respectively, yield the same 
numerical value of the contraction coefficient when the ratio 
of orifice size to tank size is the same. Thus, the two- 
dimensional analysis was extended to provide solutions for 
axisymmetric flow problems. Contraction coefficients 
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determined by Rouse and Abul-Fetouh (6) were used in the 
present investigation for analytically predicting choking 
conditions (Chapter 8) . 

A sizable compilation of many of the major references 
in two-phase flow has been given by Gouse (7) who lists 
thousands of papers in three sections for organizational 
purposes. There was no attempt by the author to discuss 
material within each individual paper, however. Rather, 
cross referencing by title, subject, and author was provided 
to facilitate the usefulness of the collection as a very 
desirable aid to research in two-phase flow. 

Among the more specific research investigations of im- 
portance was that of Frinhard and Stephenson (8) which was 
directed at temperature and scale effects of cavitation in 
two types of jets. Incipient and desinent cavitation numbers 
were presented and a correlation between superheat , -j et area, 
and flashing was given. Similar research was performed by 
Numachi, et. al. (9) who present data on effects of cavita- 
tion on the discharge coefficient for sharp-edged orifice 
plates. Kobayashi (10) performed a study similar to that 
described in Reference (9) ; however, his work on discharge 
coefficients and cavitation was performed using standard 
ASME long-radius flow nozzles and ISA flow nozzles. Gorinda 
Rao, et. al. (11) investigated cavitation in quadrant-edged 
orifices similar to those used by Numachi (9, 12), The 
researchers of Reference (11) were concerned with the 
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applications of graphical correlations between the cavita- 
tion number and coefficient of discharge for standard orifice 
meters . 

The problem of scale effects in cavitating systems has 
been described extensively by Hall and Wislicenus (13) . 

Some of the modeling parameters were tabulated by these in- 
vestigators, and their effects were illustrated by using the 
similarity principle. A significant conclusion was that 
further experimentation and analysis should be conducted. 

A similar treatment of scale effects of cavitation has been 
given by Oshima (14) . This study was restricted to the 
special case of axisymmetric bodies. Experimental data of 
Kermeen, et al. (15) and an analytical solution by Knapp 
(16) are included in Reference (14) . The main interest in 
these papers lies in the proof of validity of applying the 
cavitation number to different geometric systems. The appli- 
cation of similarity conditions to cavitation is discussed 
in References (14-16), and its validity is supported on the 
basis of experimental results. 

Hall (17) has noted a distinction between the appearance 
of cavitation under a diminishing pressure and the disappear- 
ance of cavitation under increasing pressure. Thus, he has 
defined the incipient and desinent cavitation numbers. 

Kermeen (15) and Hall and Treaster (18) came to similar con- 
clusions concerning the significance of the previous history 
of the fluid in regard to incipient and desinent cavitation. 
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Kamiyama (19) found that the longer a region of liquid was 
subjected to a low pressure, the easier became the inception 
of cavitation. 

The existence of dual cavitation numbers prompted Hall 
and Treaster (18) to consider the incipient and desinent 
cavitation numbers as being indicative of a hysteresis effect 
with a measurable cavitation-delay time. Unfortunately, the 
complexities of the problem allowed for only a qualitative, 
explanatory investigation of this effect. However, it was 
found in the present investigation that the existence of a 
measurable response time is an important consideration in the 
analytical treatment of cavitation in an oscillating pressure 
field . 

In addition to the problem of cavitation, it is appro- 
priate to consider the significance of other two-phase flow 
phenomena. Perhaps the most useful existing reference of a 
general nature is that provided by Collier and Wallis (20). 

This reference covers all aspects of two-phase flow, and 
many model formulations and specific analyses are included. 

The definitions and experimental and analytical relations 
were compiled in a systematic representation in this reference. 

A one-dimensional solution for homogeneous flow presented 
by Collier and Wallis (21) illustrates the advantages of ap- 
plying the pseudo fluid concept. Pressure drop relations and 
the one-dimensional Mach number for two -phase flow were used 
in order to analytically predict two-phase choking. The 
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general analytical formulations for two-phase floxv noted in 
Reference (21) were used in the present study in the develop- 
ment of the analyses given in Chapters 6-8. In particular, 
the influence coefficient method of evaluating and describ- 
ing the pressure gradient has been utilized in Chapter 8 to 
discuss a choking condition for two-phase flows. 

The relationship between the choking phenomenon and the 
fluid acoustic velocity for single phase flows has been well 
described by Shapiro (22) , An extension to the two-phase 
flow situation was given by Fouse (23) , whose survey reports 
surprisingly low values for two-phase acoustic velocities. 
Several analytical flow models and their associated acoustic 
velocities have been discussed by Smith (24) . Most of 
Smith's research was directed towards cryogenic fluids. 

A study of critical flow of a non-cryogenic fluid has 
been performed by Levy (25) . For the analytical model 
proposed by Levy, the velocity of each phase is represented 
by a mean value, and frictional and hydrostatic losses are 
neglected. Comparison of the predicted values with availa- 
ble measurements showed that these analytical predictions 
were in close correlation with the data. The advantage of 
this particular analysis lies in the fact that the predicted 
acoustic velocity is not dependent on the void fraction. 
Fauske (26) has used an annular type model for a critical 
flow analysis. The same model was used by Moody (27). 

Moody treated the local quality as an independent parameter 
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in his analysis. In a later publication by Moody (28), the 
existing flow models were used to predict the maximum flow 
rate for given upstream stagnation conditions. The effects 
of pipe resistance were accounted for in this analysis. An 
analytical model based upon floitf parameters in the operating 
region approaching the critical flow condition has been pre- 
sented by Nahavandi and von Mollen (29) . 

Collier and Wallis (30) list several references which 
report the use of a single phase flow equation for flow 
through a sharp edged orifice (31-34). They note that con- 
siderable uncertainty exists as to the application of the 
single phase equation when flashing occurs upstream. 

Zaloudek (33) has extended the single phase equation to 
apply to an application involving entrance choking for the 
situation in which the downstream pressure is the vapor 
pressure of the liquid corresponding to the upstream tempera- 
ture. The correlation presented by Zaloudek (33) was noted • 
as being valid for long pipes (L/d>10) with sharp edged 
entrances, as discussed by Collier and Wallis (30). A solu- 
tion which -is very similar to the choking expression of 
Zaloudek was used successfully in the present investigation 
in treating the choking of a flashing flow. The details 
of this treatment are discussed further in Chapter 8. 

Research by Fauske (33, 34) on the ’‘critical" pressure 
for short discharge pipes and subsequent modeling for this 
problem by Silver (35) were discussed further by Collier and 
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Wallis (30) . The question of whether vaporization occurs on 
the surface of the liquid core, or nucleation occurs within 
the core of liquid, is also discussed in .Reference (30). 

The existence of a critical rate implies that the flow 
achieves its sonic velocity with the possible existence of a 
shock wave. Some research in this area has been performed 
by Eddington (36) who presented descriptions for normal and 
oblique shock phenomena in a homogeneous two-phase two-com- 
ponent system at high pressures. An analysis was performed 
for isothermal conditions, and the similarity between the 
two-phase analysis and a single phase compressible flow 
analysis was discussed. Other remarks relating to the exist- 
ence of normal shocks have been made by Muir and Eichorn (37) 
and also by Elliott et al.' (38). To the best knowledge of 
the author, there has been no previous research on shock 
waves in low pressure, two-phase, single component flashing • 
flow. However, the basic approaches utilized in the analyses 
described above were most useful in the present study for 
formulating the equations governing a shock in flashing flow. 
These equations are presented in Chapter 8. 

The extension of two-phase flow research to include 
vibrational effects has been very limited in scope so far. 
Most investigations in this area have dealt with liquids 
confined in vibrating systems, and have been primarily 
directed toward the stability of the liquid-gas interface 
and the dynamic interaction between the liquid and its 
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container. Typical of work in this field is the contribution 
of Bleich (39) which consisted of an analytical solution to 
the problem of a flexible, liquid-filled tank exposed to 
longitudinal forced vibration. Boundary conditions for both 
rigid and flexible tank bottoms were considered together with 
the effects of the compressibility of the liquid. Bleich 
later extended the analysis to include prediction of the 
motion of gas bubbles in a liquid subjected to longitudinal 
vibration (40). The results have considerable significance 
with regard to propellant supply systems of liquid fuel 
rockets. Additional studies of the motions of small gas 
bubbles in vibrating liquids have been performed by Baird 
(41) and Buchanan, et al. (42). Pounder, ct al . (43) 
developed an analytical method to predict both inception and 
growth rates for gas and vapor bubbles in longitudinally 
vibrated tanks. 

Schoenhals and Overcamp (44) presented an analysis of 
the periodic pressure \ariations in a longitudinally vibrated 
liquid-filled flexible container. Their prediction of the 
oscillating pressure distribution was extended in the present 
research to yield a modified cavitation number (Chapter 7) . 

Schoenhals, Winter, and Griggs (45) have presented 
measurements which show the effect of vibration on the liquid 
discharge rates fTom tanks subjected to longitudinal vibra- 
tion. A flow retardation was found to exist due to the vibra- 
tory conditions. Pesar (46) provided further experimental 
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verification of the retardation effect first noted in 
Reference (45) . The retardation was attributed primarily 
to a lower average pressure in the discharge tube inlet* 

Research described in References (45) and (46) was 
concerned only with gravity discharge of simulated propellant 
tanks. However, these results naturally led to the expecta- 
tion that a flow retardation could also occur under vibra- 
tory conditions in the case of forced flow due to a pressure 
drop, and this topic was investigated in the present study 
(Chapters 6-8) . 

The literature discussed up to this point has dealt 
with the various phenomena as if they were separate and 
independent physical processes. It is now appropriate to 
discuss the literature describing the effects that these 
combined phenomena can produce in a total system such as a 
liquid fuel rocket, for example. The problems associated 
with a rocket system are necessarily dependent upon many 
factors which can be included in an analysis if each of 
these factors can be accurately described. 

Fashbaugh and Streeter (47) performed a transient flow 
analysis of a liquid feed system for a liquid propellant 
rocket subjected to longitudinal vibration. This paper is 
nf considerable interest in that effects of cavitation 
(onset of two -phase flow) in the pump inlet were estimated 
from experimental data. Also, the analysis was directed 
towards a specific system, and correlation between theoretical 
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predictions and experimental data are included. The longi- 
tudinal vibrations of the system were modeled as being part 
of the POGO effect. (The dynamic interaction of the rocket 
system vri'th the unsteady combustion in the motors is referred 
to as the "POGO effect".) 

An analysis by Licht and Park (48) was performed to 
determine if POGO oscillations could develop in the Saturn V 
second stage due to interaction between the propulsion system 
and the longitudinal modes of vibration. A linear mathema- 
tical model was used. The stability (lack of POGO vibrations) 
was shown to be dependent upon pump compliance, line inert- 
ance, and structural damping. The dependence of the POGO 
effect on the parameters mentioned above' was also considered 
in Reference (48) . Since the total system was treated as a 
closed loop, each segment of the loop affected the perform- 
ance of the total system. If conditions for two-phase - gene- 
ration develop in the line, the line inertance can change; 
and as a result, the overall performance of the loop can be 
affected . 

An analog simulation of the dynamic characteristics of 
a rocket propulsion system was used by Lehner (49) . In this 
work, parameters such as pump inlet cavitation compliance 
were studied so that their influence on the total system 
performance was determined, and the characteristic contri- 
butions of each specific phenomena to the instability problem 
were obtained. 
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It should be noted that the investigations described 
above were concerned with the total rocket system, and- that 
the characteristics o£ each segment in the system can affect 
the behavior of every other segment. For example, the 
extent of pump inlet cavitation is partially dependent on 
the magnitude of the combustion instabilities in the rocket 
motor. And conversely, the combustion instabilities are 
partially dependent on the extent of the pump inlet cavita- 
tion . 

The research conducted in the present investigation was 
concerned with two-phase generation and its interaction with 
a vibrating discharge line. The flow retardation effect 
described in References (45) and (46) prompted, in part, the 
initiation of the experimental program. The following 
chapter describes the experimental facility which was con- 
structed for studying this effect and for observing the 
several other phenomena previously mentioned. 
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CHAPTER 5 

EXPERIMENTAL INVESTIGATION 


5.1 Description of Flow Facility and . 

Equipment Specifications 

The experimental facility was developed with the object- 
ive of simulation of a propellant tank discharge line which 
could be subjected to a wide range of operating conditions. 
Variables such as stagnation pressure and temperature, flow 
rates, pressure losses and gradients for pure liquid flow 
and for two-phase flow were to be measured and recorded. 

Also, those parameters not dependent on the geometry had to 
be sufficiently controllable in order to systematically 
cover the desired experimental combinations of flow variables. 
The geometry of the discharge line was an independent para- 
meter. 

A closed loop fluid flow system was designed and con- 
structed. A schematic diagram of the resulting facility is 
shown in Figure 2. The test fluid is recirculated continu- 
ously through the system, thus assuring that the properties 
can be kept constant. The prime mover is a Roth two -stage 
turbine pump, Model #3751, driven by a fifteen horse-power 
U.S. Electric motor. 




Figure 2, Schematic Diagram of Closed Loop Fluid Flow Facility 
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The fluid is discharged by the pump into a. glass lined 
surge tank and then flows either through a by-pass loop or 
through the main line to the test section. The ullage 
pressure in the surge tank is variable and is controlled by 
a National dual-meter pressure regulating valve attached to 
a high pressure nitrogen tank. The pressure can be increased 
by means of the reducing valve, and the bleed valve is used 
to reduce the pressure when desired. 

The amount of flow through the test section main de- 
livery line and through the by-pass line, respectively, 
depend upon the relative settings of the nine control valves 
in the system. All the valves shown in Figure 2 are teflon 
lined ball valves manufactured by Jenkins. 

The main by-pass line consists of two inch galvanized 
steel pipe. The £low in the by-pass system can be controlled 
so that any desired portion of it is caused to pass through 
a Dupont shell and teflon tube heat exchanger. The flow 
rate through the exchanger is adjustable from zero to a maxi- 
mum of SO gpm at 150 psi. The coolant side of the heat ex- 
changer is connected to a building main water line with a 
nominal average temperature of 55° F. It acts as a variable 
heat sink to control the temperature of the test fluid. 

Once the by-pass fluid has traversed the heat exchanger 
and main by-pass, respectively, it can either flow through 
a glass lined electric heater or through another by-pass 
line around this heater. If the test loop is being operated 
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at elevated temperatures, varying amounts of the test fluid 
can be passed through the thermostatically controlled heater 
to maintain the proper temperature in the liquid in spite of 
free convection losses from the piping. Finally the by-pass 
fluid returns to the inlet of the pump. 

The test section delivery line is also connected to the 

surge tank described above. A Fischer-Porter Model 10C150S 

dual turbine flow metering system with automatic switching 

is located in a segment of straight pipe in this line. The 
* 

meter serves as a continuous monitor of the fluid flow 
through the test section. Two one foot long flexible sections 
of steam hose serve as connectors between the delivery pipe 
and the test section. The flexible connectors are required 
in order to allow the test section to be vibrated on the 
electrodynamic shaker. The fluid passing through the test 
section is returned to the manifold system upstream of the 
pump inlet where it joins the fluid flowing in the by-pass 
system. 

The shaker on which the test section is mounted is part 
of the vibrational facility shown schematically in Figure 3. 
The Calidyne Model 177A Electrodynamic Shaker is powered by 
a lv'estinghouse Industrial Audio Amplifier, Type FG-11. The 
control console consists of a combination of Calidyne Model 
231 and Model 232 systems with a Bruel and Kjaer Audio Oscil- 
lator. A feedback loop consisting of a Bruel and Kjaer type 
4334 crystal accelerometer and an Unhlotz -Dickie Model 608PS-1 
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charge amplifier for signal conditioning are used to continu- 
ously monitor the vibrational conditions of the shaker table 
and to provide for the automatic functioning of the entire 
system. A CMC Model 226A Universal Counter-Timer is used to 
measure the frequency output of the Automatic Vibration Ex- 
citer Control. The facility provides excitation at frequen- 
cies from 10 to 2000 cycles per second and delivers forces 
up to a maximum of 5000 pounds. 

Two pressure transducers are installed permanently in 
the main delivery line to the test section. The upstream 
transducer is a Pace Model P1A (0-100 psia) , and the other 
is a Pace Model P7D (0-100 psid) differential transducer. 
Signal conditioning for both transducers is accomplished 
with, a set of Carrier Demodulators Model CD10. The absolute 
transducer measures the static pressure in the delivery line. 
The upstream stagnation pressure is obtained by correcting 
for the dynamic pressure associated with the flow which is 
calculated using the measured flow rate. The differential 
transducer is used to measure the stagnation pressure loss 
across the test section. This setup is shown schematically 
in Figure 4. 

Additional instruments permanently installed in the 
flow loop are a Fischer-Porter Model 1842 Industrial Rate 
Indicator (which serves as a readout device for the turbine 
flow meters) and a Leeds and Northrup multiple range poten- 
tiometer indicator (which is used together with an 




Figure 4. Permanent Instrumentation of 



Combined Fluid Flow and Shaker Facility. 
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iron-constantan thermocouple to measure the fluid tempera- 
ture) . 

Two different test section configurations were used in 
the experimental work. Both were constructed of plexiglass 
since visual observation of the test fluid required a trans- 
parent test section both for cavitation studies and for ob- 
servations of two-phase choking. Plexiglass was chosen 
because of its transparency as well as its availability and 
relative ease of machining and fabrication. Several differ- 
ent sizes of cylindrical plexiglass tubes were used in order 
to achieve various geometries for the simulated discharge 
line. The disadvantages of this material were associated 
with its upper temperature limit of 174° F, a relatively low 
high pressure limit, and the lack of structural rigidity of 
the test section. 

The more specific measurements associated with variables 
such as oscillating pressures and pressure gradients were 
obtained with the test configuration shown in the schematic 
diagram of Figure 5. Pressure taps located along the dis- 
charge line enabled the measurement of the static pressure 
along the wall. Data for frictional losses and pressure 
gradients in flashing flow were obtained using this test 
section. A Kistler Model 701A pressure transducer was used 
for measurement of the pressure fluctuations in the fluid 
when the test section was being longitudinally vibrated. 

This was mounted in the test section with the pressure tap 
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Figure 5. Schematic Cross-Sectional View oE Specific Test 
Section . 
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Figure 5 
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located 0.125 inch from the sharp edged entrance of a 0.375 
inch diameter discharge line. The transducer signal was 
conditioned by a Kistler SN 1911 Charge Amplifier and was 
displayed on a Tektronic Type 502A dual beam oscilloscope. 

The faces of the test section were made flat and were highly 
polished to allow for undistorted flow visualization. A 
photographic study was also conducted with this test section. 

A second and more comprehensive test section, illus- 
trated in the schematic diagram of Figure 6, was used to 
investigate other geometric configurations. No additional 
instrumentation for this test section was ne.cessary since 
all data were obtained with the permanent instrumentation 
of the fluid flow loop. 

Two Sandborn Model 150 Recorders were used for the 
simultaneous recording of pressures and flow rates as 
functions of time. Each unit contained four Model 150-400 
Power Supplies and Model 150-100 AC-DC Preamplifier units 
with a Model 154-100 B Four Trace Recorder. 

Visual observation of the flow in the discharge line 
under vibrating conditions was made possible with the aid of 
a Bruel and Kjaer Stroboscope Type 4910. This stroboscope 
is syncronized with the motion of the shaker tabic by ad- 
justing the frequency signal of the Audio Oscillator installed 
in the control consol. The unit has allowance for an adjust- 
able phase angle deviation from the motion of the shaker 
table. To. the eye of the observer, the test section appears 
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Figure 6. Schematic Cross-Sectional View of Variable 
Geometry Test Section. 



stationary at any chosen position along its sinusoidal dis- 
placement path. 

The test fluid used in the flow loop was water which 
was filtered in order to remove particles having diameters 
larger than four microns. Measurements of the surface 
tension were made periodically with a Fischer Porter Surface 
Tensiomat to check the water condition in the flow loop. 

5.2 Description of Experimental Procedure 

The first step, preceeding experimental measurements, 
involved calibration of the measuring instruments. Both 
Pace transducers were calibrated by means of a dead-weight 
tester and a standard pressure gauge of sufficiently high 
accuracy. The flow meters were provided with an internally 
adjustable calibration circuit which was installed within 
the rate indicator. 

A typical liquid flow experiment was initiated by 
mounting the appropriate test section on the shaker table. 
The flow loop was then filled with water which was brought 
to the desired temperature. Finally all instrumentation 
was checked for proper operation. 

All flow control valves were initially set in the full 
open position, and the pump was started. By proper adjust- 
ment of the valves, the fractions of the total pump flow 
through the by-pass line and the test section could be pre- 
cisely regulated. Restriction of fluid flow through the 
by-pass line increased fluid flo\v r through the test section 
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and vice-versa. The pressure in the test section could be 
controlled in two different manners. The first control was 
achieved by simply changing the ambient pressure in the 
surge tank by regulating a pressure reducing valve which was 
connected between the surge tank and the high pressure nitro- 
gen bottle. In the second method, use was made of the vari- 
ous valves located in the loop. Proper adjustment of indi- 
vidual control valves permitted the selection of an appro- 
priate stagnation pressure in the main delivery line for 
any desired flow rate. With both the stagnation pressure 
and the flow rate adjustable, many different types of experi- 
ments could be carried out within the performance range of 
the pump. 

In a typical cavitation experiment, the stagnation 
pressure was held constant, and the flow velocity was in- 
creased until cavitation appeared. The stagnation pressure 
was then changed and the same operation was repeated for "this 
new stagnation pressure value. If the test section was vi- 
brated during a test run, the observer used the stroboscopic 
light to observe the cavitation, and the test section then 
appeared as a stationary object. 

A typical choking experiment was conducted in a very 
similar manner. First, choking of the flow was induced at 
a given upstream stagnation pressure by reducing the exit 
pressure to a limiting value beyond which no further increase 
in the volumetric flow rate was observed. Although a normal 
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shock could occur at various locations in the discharge line, 
the exit pressure was always adjusted to a value which would 
drive the shock into the exit plane o£ the discharge line. 

The stagnation pressure was then increased with a simultane- 
ous adjustment in flow rate to maintain the two-phase shock 
in the exit plane of the discharge line. In this manner 
many different combinations of flow rates and stagnation 
pressures were obtained. Choking experiments were performed 
in exactly the same manner when the test section vibrated. 

The following three chapters contain the results of the 
experimental work, and they also include the corresponding 
analytical predictions for most of the cases Investigated. 
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CHAPTER 6 

SINGLE PHASE LIQUID FLOW IN A DISCHARGE LINE 


6.1 Flow in the Absence of Vibration 
The classical one-dimensional relations defining liquid 
flow in a circular pipe were used to predict static pressures 
and stagnation pressure losses in the discharge line. The 
objective of this investigative phase was to obtain a compari 
son with experimental data which would serve as a criterion 
for the reliability of the facility. 

It follows from Streeter (1) that the static pressure in 
a pipe for one-dimensional flow of a single phase liquid is 
given by the relation, 



where 


^entrance ^friction ” ^dynamic + p f£ z cos ^ 

( 6 . 1 ) 

AP entrance - stagnation pressure loss 
due to vena-contracta 

APfriction = sta g nation pressure loss 

due to friction in the pipe 

^dynamic = d y namic pressure of liquid 
due to liquid motion 
P n = static pressure at location 
"n" along discharge line 
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= angle of inclination from 
the vertical 

The definitions given by Streeter (1) are listed below using 
the parameter defined as V^, the apparent liquid velocity. 
This term is used throughout the analysis, and it is careful- 
ly defined to prevent misunderstanding. The apparent liquid 
velocity is defined as the volumetric liquid flow rate mea- 
sured by the turbine flow meters in the delivery line, divid- 
ed by the cross sectional area of the discharge line. This 
means that when single phase liquid is flowing through the 
discharge line, the average liquid velocity in the discharge 
line equals the apparent liquid velocity. However, when a 
two phase (liquid and vapor) mixture is flowing in the dis- 
charge line, the average flow velocity does not equal the 

apparent liquid velocity but is a function of' the void 

* 

fraction and the term V£. Thus for the single phase liquid, 
the definitions from Streeter (1) are 


AP 


entrance 


“ k (t Of v f ) 


( 6 . 2 ) 


where K * 0.5 for the geometry investigated. 

z 


AP 


friction 


* 2 , 


- r £ p £ V D 


(6.3) 


1 * 2 
^dynamic 2 p f 


(6.4) 


Substitution of the definitions into equation (6.1) results 
in the one dimensional equation governing the static pressure 
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in the discharge line. 


P n = p s + PfS z cos ' 


£z 

d” 


+ 1 


51 * 
31 2 


P £ 


V, 


*2 


(6.5) 


Friction was considered to be important not only for the 
single phase analysis but also for the two-phase analysis to 
be discussed later. Shames (2) presents friction factors in 
a Moody diagram. Representative values for a very smooth 
pipe are listed in Table 1. 


Table 1 

Friction Factors 


(ft/sec) 

Reynolds No. 

f 

10 

33,800 

0.0222 

20 

67,600 

0.0195 

30 

101,400 

0.0178 

40 

135,000 

0.0170 

50 

109,000 

0.0163 

60 

202,500 

0.0159 

70 

236,000 

0.0154 


An experimental verification of these values was per- 
formed. Tv/o different pressure transducers were connected 
to pressure taps numbers 2 and 3, and numbers 3 and 4 respect- 
ively as shown in the schematic diagram, Figure 7. The 
pressure difference between the stations under flow conditions 
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Figure 7. System Configuration of "Specific" Discharge Line 
and Propellant Tank, 
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represented the frictional pressure loss. A comparison 
between experimental data and predictions generated with 
equation (6.3) using friction factors from Table 1 is shown 
in Figure 8. The results of the one dimensional equation 
yielding static pressures in the discharge line, equation 
(6.5), are compared with experimental data in Figure 9, 
where the static pressures p n , n - 1, 2, 3, 4, (the subscripts 
1, 2, etc. define locations along the discharge line) are as 
illustrated in Figure 7. The value for the Reynolds number 
listed in Table 1 was obtained by using property values for 
water at 80° F. The majority of experiments for this initial 
phase of the investigation were conducted at this temperature. 

An important feature of the geometric configuration 
illustrated in Figure 7 is the stagnation pressure loss 
caused by the sudden expansion of the single phase fluid as 
it discharges from the simulated discharge line. Ihe flow 
pattern during expansion is shown schematically in Figure 10 
which was taken from Kenyon (3) . Streeter (1) derives an 
expression for the stagnation pressure loss due to a sudden 
expansion for the stagnation pressure loss due to a sudden 
expansion in terms of characteristic diameters which are 
indicated in Figure 10. The expression for the expansion 
pressure loss is noted subsequently. 

Ap . = stagnation pressure loss 

expansion 

due to sudden expansion 

AP - = (1 - d 2 / D 2 ) 2 i p f V*2 

u expansion ^ ' L * 1 


( 6 . 6 ) 



FRICTIONAL HEAD LOSS 2 Ap 24 = p 2 - p 4 (PSI) 



FRICTION FACTORS LISTED 
IN TABLE I 

EQUATION (6.3) _ 


P s SYM 
40 ~o~ 
60 □ 
80 A 


20 40 

V f * (Ft/sec.) 


Figure S, Illustration of Measured Frictional Head Loss 
as a Function of Apparent Liquid Velocity. 





Figure 9. Illustration of Static Pressure Distribution Along Discharge Line Wall 
with Different Apparent Liquid Velocities at One Constant Stagnation 
Pressure. 



D 


SUDDEN CONTRACTION 



SUDDEN EXPANSION 


Figure 10. Single Phase Liquid Flow Through Sudden Area 
Changes from Kenyon (3) . 
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The total sum of all stagnation pressure loss terms 
represents the irrecoverable pressure lost by the fluid 
flowing through the discharge line. 


AP 


loss 


AP 


loss 


^entrance + ^friction + ^expansion 
[0.5 + + (1 - dVD’p] | 0f vj* 05.7) 


This pressure loss was measured by the differential 
pressure transducer connected to the pressure taps located 
between the test section delivery line and exit line which 
are shown in Figure 4. The signal generated by the trans- 
ducer was corrected for the losses caused by the elbows and 
the friction in the test section by measuring these losses 
without a discharge line present in the system. A comparison 
of experimental data for the stagnation pressure loss across 
the discharge line with the results of equation (6.7) is pre- 
sented in Figure 11 for several different geometries. 

The expansion pressure loss relation is assumed to be 
valid for the discharge line configuration shown in Figure 6 
as well as for that configuration shown in Figure 5. This 
assumption is valid because of the nearly identical values of 
the loss coefficient, 

(l - d 2 /D 2 ) 2 * 1.0 (6.8) 


The condition of d<<0 permitted the flow to be treated as 
that of a free jet expansion problem. 
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Figure 11. Stagnation Pressure Loss Versus Apparent Liquid 
Velocity for Different Liquid Flow Geometries. 
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The one-dimensional analysis for the static pressure 
distribution presented above becomes inadequate in the 
proximity of the sharp edged corner of the entrance to the 
discharge line. The experimental values of the static 
pressure, pj , shown in Figure 9 were found to be very dif- 
ferent from the values predicted by equation (6.5). The 
reason for this discrepancy was obviously the vena-contracta 
effect illustrated in Figure 10. In order to predict the 
static pressure, p A , at the entrance to the discharge line 
and the local static pressure across the entrance plane, a 
three-dimensional analysis would have to be performed. 

Order of magnitude calculations demonstrated that the pre- 
dictions obtained for a one-dimensional model differed by 
an order of magnitude from the experimentally determined 
pressures . 

Since no three-dimensional solution to this problem was 
found in the literature and since a solution of this nature 
was beyond the scope of this investigation, it was decided 
to employ a less sophisticated two-dimensional model. The 
method of conformal transformation could have been applied 
to this type of two-dimensional flow problem. However, this 
method of solution was not considered. A conformal transfor- 
mation flow solution to the geometry of the discharge line 
necessarily required negative pressure values at the sharp 
edged entrance. It also required that the fluid did not 
separate from the physical boundary. The vena-contracta 
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effect combined with the physical reality of having a finite 
pressure at the corner were powerful arguments against using 
a boundary capping solution for the velocity and pressure fields. 

A free-streamline model was considered as being more 
representative of the fluid flow near the corner. Although 
the free streamline concept is usually applied to a liquid 
in contact with a fluid of much lower density, nevertheless 
the free streamline was considered adequate for the descrip- 
tion of the flow in the region at the entrance to the discharge 

line. This assumption was based on two arguments: 

z 

1. The area ratio (d/D) was always very small, and 
the flow rates were usually high enough to cause 
significant acceleration forces on the fluid 
flowing into the discharge line. 

2. The region of application was very limited in 
size . 

Since the theory of the free-streamline can be found in 
many textbooks, it is not rederived here. For a general 
solution to the problem of flow from a vessel, von Mises (4) 
showed that most flow situations were particular cases of a 
well defined problem. This problem is illustrated in Figure 
12 which was taken from von Mises (4) . For the problem 
under investigation, shown in Figure 13, the constants used 
by von Mises are given below. 

The values of von Mises’ constants are: 

K « 1/2 . 



e = o 

a = b 
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j = D/2 

A and E coincide 

The free streamline analysis lead to the determination 
of the contraction coefficient which is defined by the 
asymptotic width of the resulting jet. 

C c = contraction coefficient 

C = asymptotic area of free jet 

u c “ area of orifice 

C c = (3/b in Figure 13 (6.9) 

Since the velocity along the free streamline is constant 
and the final flow area is known, the pressure along the free 
streamline can be determined immediately from Bernoulli's 
equation. 

Table 2 represents the tabulations of contraction coef- 
ficients arrived at by the solutions of the two-dimensional 
problem posed in Figure 13 and cited by Gurevich (5) . 


Table 2 

Contraction Coefficients 


b/B 

c G 

b/B 

C c 

O 

o 

0.611 

0.5 

0.644 

0.1 

0.612 

0.6 

0.622 

0.2 

0.616 

0.7 

0.687 

0.3 

0.622 

0,8 

0.722 

0.4 

0.633 

0.9 

0.781 


» j 


t PLANE 




Figure 12. General Free Streamline Liquid Flow Problem and 
Method of Solution from von Mises (4) . 
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Figure 13. Two-Dimensional Liquid Flow Problem from 
von Mises (4) . 



Figure 14. Axisymmetric Liquid Flow Problem from Rouse and 
Abul-Fetouh (6) , 
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The tabulations of von Mises are solutions for two- 
dimensional problems only. The problem on hand is a special 
case of a three-dimensional flow, namely axisymmetric flow. 
An axisymmetric flow problem is formulated helots. 


!i „ 1 « v 

9z r 9r z 


( 6 . 10 ) 


9<t> _ _ 1_ 9>| ) _ y 

3r r 3z " r 


( 6 . 11 ) 


where $ = potential function 

= stream function 
r = radial coordinate 
z = longitudinal coordinate 
0 = angular position coordinate 
By differentiating and combining equations (6.10) and (6.11), 
Laplace's equations for axisymmetric flow are obtained 

+ lit + 1M = o (6.12) 

al 7 ar 2 r 3r 

i + ili - I 14. » o (6.13) 

3z 2 ar 2 r 3 r 

The hydrodynamic problem is solved when the functional 
relationship of either <}>(r,z) or ^(r,z) is determined. Along 
the free streamline, the boundary conditions are' known. 


= 0 (6.14) 

an 

ip = constant 


where n = outward drawn normal 
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The solution of the above equations leads to great 
mathematical difficulties. While the only known three- 
dimensional solutions are given for axisymmetric jet flows, 
Gurevich (5) points out that papers dealing with axisymmetric 
jet flows either describe theorems of a general or qualita- 
tive nature or present solutions which are at best numerical 
or approximate. 

An analogue solution of the problem illustrated in 
Figure 14 was given by Rouse and Abul-Fetouh (6). These 
authors used an electrohydrodynamic analogy to solve the 
problem of flow of a circular jet from an infinitely long 
circular vessel. Rouse and Abul-Fetouh (6) have tabulated 
the contraction coefficients of the axisymmetric flow problem 
which was solved earlier ( for a two-dimensional model by 
von Mises (4). The coefficients are listed in Table 3. 


Table 3 

Comparison of Contraction Coefficients for Slots and Orifices. 
(Taken from Rouse and Abul-Fetouh (6) and von Mises (4).} 


(b/B) 2 = (d/D) 2 

C c , Slot 
von Mises 

C c , Orifice 
Rouse and Abul-Fetouh 

0.00 

0.611 

0,612 

0.25 

0.644 

0.644 

0.50 

0.689 

0.691 

0.75 

0.757 

0.757 

1.00 

1.000 

1.000 
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At equal values of the ratios describing the geometries 
of the problems, 



the contraction coefficients are equal. 

Using the theory of free streamline flow and the fact 
that the contraction coefficient yields an effective flow 
area, the static pressure in the area between the free stream- 
line and the wall of the discharge line geometry illustrated 
in Figure 10 can be approximated. Since the pressure and 
velocity along this particular streamline are constant from 
its starting point at the sharp-edged entrance to its final 
point of deterioration, Bernoulli's equation can be used to 
obtain the value of the static pressure along this streamline. 

Neglecting the gravity force for the liquid within the 
small entrance region, the total pressure along the stream- 
line is defined as: 

p s - p’ + j P f V' 2 (6.17) 


where p' h static pressure along the free streamline 

V' = constant velocity along the free streamline 
Using the contraction coefficient and the continuity 


equationtwith equation 

(6.17) the static pressure 

equation was 

obtained as follows. 



Continuity : 

Pl A l V l “ P 2 A 2 V 2 

(6.18) 

»£ V £* A d 

= ^ V * A iet = Pf V ' C c A d 

(6.19) 


or V’ = V-*/C„ 
f c 

(6.20) 



STATIC PRESSURE , p, (PSiA) 


P s =80 


l-D BERNOULLI’S EQUATION 
PREDICTIONS , 


P s =60 


P s = 40 


-IPs 

S' 

y~ 

go 

40 

o 

- 60 

€> 

80 

O 


EQUATION (6.21) 
PREDICTIONS ^ 


k v 
V f (Ft. /sec.) 

Figure 15. Static Pressure as a Function n-F a 

S lo fty 1/3 Diameter DownstJeam ij quid 

. Edged Entrance _yw«stream ±rom the Share 

Stagnation Pressures? DlSChar ^ e iine Under Various 
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V f * (Ft. /sec.) 


0 20 40 60 80 



Figure 16. Dynamic Pressure as a Function of Apparent Liquid 
Velocity 1/3 Diameter Downstream from the Sharp 
Edged Entrance in the Discharge Line Under Various 
Stagnation Pressures. 
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Therefore 

p' - p s - I Pf( v f‘ /C c ) 2 (6 - 21 ’ 

.Equation (6.21) describes the static pressure along the 
free streamline only as long as the jet persists. The region 
noted as "eddies” in Figure 10 represents the limit of the 
applicability of equation (6.21). The static pressure in the 
fluid in the eddy flow region is very difficult to predict 
analytically. 

Data taken by a static pressure probe located 0.125 inch 
downstream from the entrance of a 0.375 inch diameter dis- 
charge line are compared in Figures 15 and 16 with the analy- 
tical predictions obtained by solving equation (6.21). Note 
that in these figures, the relation (6.21) is compared to a 
one-dimensional Bernoulli's equation. The only factor con- 
sidered in the one-dimensional equation was the- area change 
in going from the flow area with diameter, D, to the discharge 
line with diameter, d. 

6.2 Single Phase Flow in the Presence 
of Longitudinal Vibration 

In order to study the influence of vibration on the 
single phase flow in a discharge line, an experimental in- 
vestigation was conducted. The influence of frequency and 
acceleration level was measured for various flow rates and 
geometric flow configurations of the discharge line. While 
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vibrations in a real situation are usually random in nature, 
the correlation between frequency and acceleration and 
physical flow effects could best be determined for a sinu- 
soidal vibration situation. This restriction to sinusoidal 
vibrations was also necessitated by the lack of random vibra- 
tion capabilities of the vibrator facility. 

Initial data were taken for the pressure loss across 
the test section with the discharge line removed. In this 
calibration procedure the pressure losses caused by the 
elbows and other fittings in the test section were determined. 
Thus data obtained with the discharge line installed in the 
test section could be corrected for the pressure losses not 
directly associated with the discharge line. 

Data obtained for the stagnation pressure loss across 
the test section in the absence of vibration with the dis- 
charge line removed are presented in Figure 17. A comparison 
for cases where frequency and acceleration are independent 
variables is illustrated in Figure 18. Data are shown and a 
comparison with the best fit curve for the non-vibratory 
situation (F - 0) is included in these plots. 

The obvious conclusion is that the stagnation pressure 
loss increases where vibration is present. For a given 
frequency, the pressure loss increases with increasing ac- 
celeration level. This is consistent with the fact that at 
higher acceleration levels for a given frequency the dis- 
placement amplitude is larger than at lower acceleration 
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levels. Thus the larger magnitude of vibration would cause 
a larger loss in stagnation pressure. 

The results obtained in the experimental investigation 
of the stagnation pressure losses across the test section 
with the discharge line installed are similar to those 
described above. A vibrated test section where the discharge 
line was installed showed an increase in the stagnation 
pressure loss. However a significant difference was noted 
when the data was corrected to represent the stagnation 
pressure losses only across the discharge line rather than 
-the combination of the line plus test section. When the 
values for the pressure losses across the test section 
without discharge line present were subtracted from the 
pressure losses across the test section with the discharge 
line installed, the pressure loss values were lower for the 
cases where vibration was present than for the non-vibrating 
case. This can be stated in a different manner by saying 
that the pressure loss increase for the test section without 
the discharge line was greater than the pressure loss increase 
for the test section with the discharge line installed under 
the same conditions of flow rate, vibrational frequency, and 
acceleration level. These data are presented in Figures 19 
and 20 for various vibration conditions and geometries. 

The fact that vibration of a discharge line retards the 
flow has been noted previously by Pesar (46) and in reference 
(45) . This corresponds to having a larger pressure loss at 




Q (PER CENT OF MAXIMUM) 

Figure 17, Stagnation Pressure Loss Across the Test Section in Absence of Discharge 
Line Given as a Function of the Volumetric Flow Rate, Q(Qmax = 415 GPM) . 
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Figure 18. 


Stagnation Pressure Loss Across the Test Section 
in Absence of Discharge Line Given as a Function 
of the Volumetric Flow Rate, QCQmax * 45 GPM) , for 
Various Vibrational Conditions. 
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. Stagnation Pressure Loss Across the Discharge Line 
for Various Vibratory and Non-Vibratory Conditions 
as a Function of Apparent Liquid Velocity 
(d * 0.375 in. , L = 24 in.) . 


Figure 19 
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Figure 20 . 



Stagnation Pressure Loss Across the Discharge Line 
for Various Vibratory and Non-Vibratory Conditions 
as a Function of Apparent Liquid Velocity 
(d = 0.37S in., L = 12 in.). 
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a given flow rate which indeed the data presented in Figure 
18 demonstrate. 


6.3 Conclusions 

At this phase of the investigation, the check-out of 
the experimental facility was completed. The experimental 
facility was proved capable of providing those fluid flow 
conditions necessary to investigate high speed liquid flow 
in a simulated discharge line. Comparison of experimentally 
obtained data with the well known predictions for single 
phase fluid flow showed excellent agreement. The one-dimen- 
sional flow equation, (6.21), for the stagnation pressure 
loss across the discharge line evidently described the flow 
to a high degree of accuracy as exhibited by the comparisons 
with the experimental data. 

The objective of finding a model describing accurately 
the multi-dimensional flow pattern in the entrance of the 
discharge line was achieved by the application of the free 
streamline theory. The theory of free jets combined with the 
axisymmetric solution presented by Rouse and Abul-Fetouh (6) 
enabled a one-dimensional analytical model to be developed 
which was verified by experimental results. 

The experimental investigation of single phase fluid 
flow subjected to the influence of longitudinal vibration 
showed that stagnation pressure losses increased with increas 
ing acceleration levels. The frequency influence was less 
pronounced; however, in the high frequency range, the 
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pressure losses tended to decrease in magnitude. An analy- 
tical solution' to this problem was not attempted. 

The inconsistency of having smaller stagnation pressure 
losses across the discharge line can be attributed, perhaps, 
to difficulties in proper instrumentation and signal condi- 
tioning in the pressure transducers. Resonances within the 
tubing leading to the sensing elements inside the transducers 
could have produced errors not readily noticeable or easily 
eliminated . 

Both the one-dimensional and the multi-dimensional solu- 
tions derived for single phase flow in discharge lines were 
now extended into two-phase flow situations. This is dis- 
cussed in the next chapter. 
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CHAPTER 7 

CAVITATION INCEPTION IN A DISCHARGE LINE 


7.1 The Flow Problem in the Non-Vibrating Case 
Cavitation is defined as the formation of vapor filled 
voids in a liquid caused by a local pressure reduction as a 
result of the dynamics of the flow situation. Potential 
flow theory for a pure liquid predicts that a minimum pressure 
will occur on a solid boundary. Cavitation in a pure liquid 
will occur at some given value of the ratio. 


a - 


P 


0 " P 

1 

2 p f 


minimum 

V 


(7.1) 


where p 0 = reference static pressure 

and V 0 = a reference velocity usually taken at the same 
point as p 0 

For a liquid which generally can not support a tension, 
cavitation will take place when the minimum static pressure 
equals the vapor pressure. Thus the classical definition of 
the cavitation parameter can be formulated. 


a - 


Po - Py 

I ” TJ Z 

7 p f V o 


(7.2) 


Similarity requires that the cavitation numbers in similar 
systems be 'identical at thd point of incipient cavitatioii'. 
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In a real liquid, consideration of the surface tension 
of the liquid leads to a modification of the assumption that 
cavitation appears as soon as the local pressure reaches the 
vapor pressure of the liquid. It should be noted that for 
a single spherical bubble (i.e. cavity), the concept of 
thermodynamic equilibrium has to be used in order to estab- 
lish a criterion for the existence of a cavity. A balance 
of forces has to take into account the surface tension of 
the liquid. This balance leads to the well known relation- 
ship , 

” r o 2 (P v - P*} = 2*r 0 Y (7.3) 


where r 0 = bubble radius 

Y - surface tension 

The resulting pressure difference between the liquid and the 
vapor is given by'. 


P 


v 



(7.4) 


Inspection of equation (7.4) shows that the pressure inside 
the bubble must be larger than that of the liquid for a 
bubble to exist. Since a basic assumption is that of thermal 
equilibrium, the liquid must be superheated to allow the 
vapor bubble to exist. 

It follows from equation (7.4) that as the radius of 


curvature of the bubble becomes larger, the amount of super- 
heat required for equilibrium becomes smaller. For a situa- 
tion where the radius of curvature of the liquid vapor 
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interface could be considered very large, the pressure dif- 
ference and amount of superheat both go to zero. 

Application of the definition of the cavitation number 
as given in equation (7.2) leads to the conclusion that 
cavitation would occur under a well definable flow situation, 
provided this flow situation can be treated adequately with 
an analytical model. As discussed earlier, this prediction 
is quite difficult in the entrance region of the discharge 
line under investigation. If the flow velocity through the 
discharge line is increased while a constant upstream stag- 
nation pressure is maintained, then cavitation will even- 
tually occur when the static pressure along the free stream- 
line reaches the vapor pressure of the liquid. Initially 
this cavitation begins at the sharp edge of the discharge 
line entrance. The analysis presented in Chapter 6 gives 
the value of the free streamline pressure as a function of 
the velocity, stagnation pressure, the geometric parameter 
and its resulting contraction coefficient. 

Using the definition of the cavitation number given by 
equation (7.2), a modified cavitation number which is inde- 
pendent of the geometry of the discharge line was derived 
and used for this special type of problem. The modified 
cavitation number is being introduced in order to express 
the cavitation parameter in terms of the apparent liquid 
velocity, V^, which proved quite useful earlier in the 
analysis and experiments. The continuity equation can be 
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expressed in terms of the. flow at the reference section 
(i.e. the location where p 0 and V 0 are measured in the 
delivery line) and the flow inside the discharge line assum- 
ing a uniform velocity profile. 

Continuity : 

p £ A o V o = p f A d V f* ( 7 *. 5 ) 

V 0 = — V * (7.6) 

0 V 

Substituting equation (7.6) into equation (7.2) leads to 


Po - P y 

7 Pf v / 2 


(7.7) 


Now the definition for the modified cavitation number, 
is stated. 


a* = 



a 


a* 




(7.8) 


Substituting the cavitation number, a, as formulated in 
equation (7.7) into the above definition yields a simple 
relation for the modified cavitation number, a*. 


a* = 


Po 

1 

2 p f 



(7.9) 


The values of a* are obtained experimentally from the 
static pressures, p 0 , which are measured in the delivery 
line, and from the value of which is obtained from the 


turbine flow meter. 
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The theoretical value for a* can be derived from the 
free streamline analysis presented in Chapter 6. Bernoulli's 
equation along a streamline at the reference location is 
written in terms of the stagnation pressure 



(7.10) 


The pressure along the free streamline was defined in equa- 
tion (6.21) and is repeated here. 


Pi “ P s " 2 p f (V /C c) 


( 6 . 21 ) 


Given the constraints that the pressure along the free stream- 
line equals the vapor pressure and cavitation occurs, then 
equating of equations (6.21) and (7.10) leads to the relation, 


1 . r 2 

2 p f V o 


Substituting equation (7.6) into the above equation yields 



(7.12) 


Now substitution -of equation (7.12) into equation (7.9) leads 
to the final form of the modified cavitation number, a*. 


. (7.13) 


Equation 7.13 is an exact relation. For systems where 
D Q >>d, the approximation for the area ratio noted below is 


valid. 
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It becomes: 

— << 1 (7.14) 

V 

The constraint exemplified by equation (7.14) leads to a 
simplification of the modified cavitation number. 


a* = 


V 


(7. IS) 


The validity of the approximation for the modified 
cavitation number shown in its final form by equation- (7 . 15) 
is illustrated by using the contraction coefficient listed 
by Rouse (6) for the ratios of diameters used in the work. 

For the following conditions, 

§ < O.i 


C c = 0.611 

the corresponding value for the modified cavitation number 
is 

a* * 3.72 (7.16) 


Since the contribution of the fourth power of the ratio d/D 0 
is less than 2 % over the entire range of delivery line geo- 
metries used in the investigation, its contribution to the 
modified cavitation number can be neglected. Thus the modi- 
fied cavitation number becomes. 


o* = 


fd } 


C c a 


D 


o' 


C c 2 


(7.17) 


Typical experimental data to compare with the theoretical 


prediction of (7.17) were obtained by keeping the upstream 




V* (Ft. /sec.) 


Figure 21. Modified Cavitation Number as a Function of 

Apparent Liquid Velocity for Various Discharge 
Line Geometries. 
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stagnation pressure constant and increasing the flow velocity 
until cavitation appeared. A more complete description of 
the experimental procedure used’ in this phase of the inves- 
tigation is presented in Chapter 5. 

A comparison between experimental data and a value 
calculated with equation (7,17) is illustrated for several 
geometries in Figure 21. 


7.2 Cavitation Inception in a Discharge Line under 
the Influence of Longitudinal Vibration 
It has already been shown that the inception of cavita- 
tion is a function of the local static pressure. An inves- 
tigation by Schoenhals (44) of a flexible liquid-filled 
cylinder under the influence of longitudinal vibration proved 
that the local pressure oscillation could- be predicted 
reasonably well. A schematic diagram of the system investi- 
gated in reference (44) together with a schematic diagram 
of the test section used in this investigation are shown in 
Figure 22. Referring to the schematic diagram for definitions 
an expression for the first natural frequency of the combined 
tank and liquid system was derived and presented by Schoenhals 
(44) to be 


0 ) * l 


IT 

2 


(7.18) 


where C = acoustic velocity of the system. 

Although Schoenhals (44) derived the above expression for a 
liquid filled closed container, it was applied successfully 
to the test section used in this investigation. This proced- 
ure was- justified in the light of the results obtained in a 
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concurrent investigation by Pesar (46) , who found that for 
small flow velocities, the natural frequency' for a given 
system is independent of the flow velocity. The acoustic 
velocity, C, is dependent upon the structural parameters of 
the tank, the acoustic velocity in the water, , and the 
bulk modulus of plexiglass of which the tank was fabricated. 
This derivation is presented in total in the reference; the 
final form of the expression is given below. 


C 


1 



TJT 


(7.19) 


Using equation (7.18), a resonant frequency of approxi- 
mately 85 cps was calculated for the test section depicted 
in Figure 5. The length, l, in the equation represents the 
distance from the top of the elbow on the top of the test 
section to the entrance plane of the discharge line. Under 
conditions of no flow, a comparison of analytical predictions 
given by Schoenhals (44) and experimental data obtained with 
the test section was made. 

Pressure measurements were obtained with a Kistler trans 
ducer as indicated in the schematic diagram of Figure 5. A com 
parison of the experimental data obtained in this phase of 
the investigation with the theoretical prediction of pressure 
amplitudes as a function of frequency is shown in Figure 23. 
The data are somewhat scattered and are lower than the calcu- 
lated values of the oscillating pressure amplitudes. However 
the prediction does give the trend and order of magnitude of 
the data. 
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Figure 22. Schematic Diagram of Liquid Filled, Flexible 
Cylinder from Schoenhals (44) Compared to a 
Schematic Diagram of the Test Section Used in 
This Investigation. 

In a subsequent investigation, data of pressure ampli- 
tudes were obtained for the same values of the stagnation 
pressure but now under moderate flow conditions, i.e. moderate 
apparent liquid flow velocities. Again the experimentally ob- 
tained pressure amplitudes as a function of frequency under 
various stagnation pressure conditions are compared with 
Schoenhals' theoretical prediction. Fairly good agreement of 
the experimental data with the theory,, an observation which 
was also made by Pesar (46), is shown in Figure 24. 
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Under vibratory conditions, the stagnation pressure in 
the discharge line is composed of a steady state value upon 
which a pressure fluctuation is superimposed caused by the 
oscillating motion of the system. The instantaneous value 
for the stagnation pressure is given by, 

p s (t) = P g + AP g cos to t (7.20) 

where AP g is a known function of frequency and acceleration 
level as given by Schoenhals (44) . No phase angle relation- 
ship is shown in equation (7,20) since the discussion is 
concerned with the value of P (t) at the entrance plane of 
the discharge line where the phase angle value would be zero 
if included. By using the expression for the static pressure, 
p’ , discussed in Chapter 6 as equation (6.21), the fluctuation 
Ap' of the static pressure along the free streamline is shown 
to be equal to A p s . 

Neglecting the hydrostatic pressure term, Bernoulli's 
equation is valid along the free streamline in the -entrance 
plane of the discharge line. 


P s = P' + 2 Pf 


VI 2 


(7.21) 


This equation can be modified to account for the periodically 
oscillating pressure. 

fv^r 

(7.22) 


P s + AP S cos wt = p'(t) + ^ Pf 


The above equation follows from the approximation that only 
the static pressure is subjected to periodic fluctuations. 



Figure 23. 



F (CPS) 

Pressure Oscillations in the Discharge Line Under No Flow Conditions 
Compared to the Prediction by Schoenhals (44) . 








F(CPS) 

Figure 24. Pressure Oscillations in the Discharge Line Under Flow Conditions 
Compared to the Prediction by Schoenhals (44) . 
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It will be shown that the periodic velocity fluctuations 
which are induced by the vibration are negligibly small. It 
follows from the base motion of the test section shown in 
Figure 22 that displacement, velocity, and acceleration of 
the system under consideration obey the relations noted as 
equations (7 . 23) . 

... v iwt 
y(t) = Y l e 

y (t) = Y L ia)e lwt = ive 1U)t (7.23) 

y (t) = - Y L w 2 e iwt . -ae iwt 

Considering the magnitudes of the velocity and acceleration 
expressions, the magnitude of the oscillating velocity is 
expressed as a function of acceleration and frequency. 

i v 1 = | a/w | (7.24) 


The maximum value for the velocity magnitude at a given 
acceleration level occurs at the lowest frequency. For this 
investigation, this minimum frequency was 10 cps (20n .rad/sec) 


v lma X " l a/ “ninl " l a/20l, l 


(7.25) 


A numerical evaluation of equation (7.25) illustrates why 
this oscillating velocity component can be neglected when 
compared with the minimum value of V^* . For example 


* 

f 

v| 

* 

f 


max 

min 


30 ft/sec and G 
= | 32 . 2 / 2 0 it j = 

= 30 ft/sec >> 


= 1 |a| - 32.2 ft/sec 

0.512 ft/sec 


(7.26) 


= 0.512 ft/sec 
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With the above constraints, equation (7.22) implies the 
validity of the equation 

P’ - P’(t) ~ p' + AP g cos tut (7.27) 

Equation (7.27) can be used as an argument that the definition 
for the cavitation number in a non-vibrating flow can be used 
for this particular problem upon one modification. That is, 
the resulting expression must include a time dependent term 
accounting for the time dependency of the static pressure 
along the stream line under consideration. The second modi- 
fied cavitation number, o*(t), is therefore the time dependent 
function 


c* ( t ) - -i- 

C 2 
c 



+ 


AP S cos tot 

1 n V * 

7 p f f . 


(7.28) 


Ihis expression can be rewritten by defining a Act* term. 


Act - 


AP 

s 

\ Pf V 


(7.29) 


Therefore the final form of the time dependent 'modified 
cavitation number is written as 


D 4 
0 


o* (t) - JL 

C 2 


+ Act* cos tot 


(7.30) 
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The term Ao* needs to be discussed firstly in its implica- 
tions on a hypothetical fluid and secondly for its implica- 
tions on a real fluid. 


Consider a hypothetical fluid which is not capable of 
sustaining a tension. Then for such a fluid it can be stated 
that the static pressure to which it might be subjected due 
to the dynamics of the flow process has to obey the relation 

’ T) . > P 

K static — vapor 

The static pressure in the fluid was measured at station 
1 as indicated in Figure 6 and its measurements were displayed 
as a function of time. For given experimental conditions, 

P s = constant 
F = constant 
G = constant 
AP g = constant 

the flow rate was increased by lowering the exit pressure in 
the test section. The static pressures, p x , and the apparent 
liquid velocity, V£* , are shown in Figure 25 as a function of 
time with the assumption that the fluid is hypothetical. A 
complete description of Figure 25 is given below: 

Case CD V f *(i) = 0 

Pi(t) « P S (D 

Under no flow conditions, the stagnation pressures 
and static pressures are identical. 
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Case (ii) V^*(ii) > V.£*(i) = 0 

p^mean) < P s (mean) 

PiCt) > P v 

The flow velocity causes a reduction in the mean 
static pressure. The reduction in the mean static 
pressure is not largo. 

Case (iii) V^*(iii) > V£*(ii) > 0 

P y < PjCmean) < P s (mean) 

(p j (mean) - P y ) < AP g 

The flow velocity has increased to a value such 
that cavitation would occur in part of the sinu- 
soidal pressure cycle if the liquid could not 
sustain any tension whatsoever. 

Case (iv) V£*(iv) > V£*(iii) >> 0 

p , (mean) t heor . < P v 
AP s ’ ( p v ' P , Cmean] ) 

The flow velocity has increased to a value such 
that the theoretical mean static pressure is below 
the vapor pressure. This case was not observed 
experimentally. Cavitation occurs almost through- 
out the cycle except foT a very small increment 
of time. 

In order to discuss the significance of Figure 25, it 
is necessary to define an additional parameter which 



(i) P, 


tor 0 


p (MEAN) = P< 

I 


(ii) 


V* > V ^ 

l v f(iir v f(i) 


P v <p (MEANXPg 


|V f(iii)> V f(ii)| 


| P v < p. (MEAN) < P s 


(iv) 


p. (MEANXP V 


! IV f*iv) >Vf1iii) 


PD = 180 270 0 90 t 

Figure 25. Pressure Oscillations in a Hypothetical Fluid as 
a Function of Time and Apparent Liquid Velocity. 
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references the point along the sinusoidal path of the shaker 

motion which is being considered. This reference parameter 

is called the phase deviation, and it is defined below. 

PD <= phase angle relationship of the 
shaker table when the strobe would flash. 

PD = 0 is the lowermost position of the 
stroke. 

PD = 90 is the neutral point on the path 
of the table with an upward velocity component 

PD = 180 is the uppermost position of the 
stroke 

PD = 270 is the neutral point on the path 
of the table with a downward velocity component 

Utilization of this parameter, PD, is made when considera- 
tion is given to the pressure fluctuations noted in Figure 25. 

The fact that the cavitation number obtained when viewing 
Case (iii) at PD =l&0 is different than the case of no vibra- 
tion is obvious as the no vibration case would need a higher 
V before cavitation would occur. Similarly but in an op- 
posite sense an observer would not detect cavitation in the 
oscillating test section if he were to observe the discharge 
line under Case (iv) at PD = 0 ; however, if the test section 

were not vibrating, cavitation would definitely be occurring. 

At this point, reference should be made to the experi- 
mental procedure adhered to in the measurement of the incep- 
tion of cavitation. The onset of cavitation was observed 
visually. A Breul and Kjaer Stroboscope Type 4910 with 
control unit was used for the investigation of cavitation. 

The features of the Stroboscope are too numerous to be 
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described here; however, one feature of which extreme use was 
made has to be mentioned. This was the phase deviation 
control which permitted synchronization of the flash from the 
lamp with the motion of the shaker table. In this manner, 
the shaker table and the test section mounted on it were 
visually static to the eye of the observer. 

In these experiments, extensive use was made of an ob- 
server to detect the onset of cavitation in the discharge 
line. The operator manipulated the controls of the fluid 
flow loop and of the electrodynamic shaker while the observer 
concentrated his attention to the cavitation in the line. 

Data could be taken for only brief periods because of the 
strain on the eyes of the observer. 

The data obtained for different PD values and presented 
in Figures 26 through 33 do not display a significant change 
in the value of o*(t). The reference static pressure, p 0 , for 
all conditions in these figures was kept at a constant value 
for each datum and ranged between 30 and 80 psia. 

Cavitation measurements were performed in both test 
sections. Since the natural frequencies of the two test 
sections were expected to be different, an additional cal- 
culation of the magnitude of the pressure fluctuation was 
made for the second test section. The natural frequency 
for the test section shown schematically in Figure 6 was 
calculated to be 210 cps . 
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As described earlier, for all values of frequency well 
below this natural frequency, i.e. F<150 cps , the magnitude 
of the total pressure fluctuation should be relatively inde- 
pendent of the frequency and only a function of the accelera- 
tion level. 

In order to obtain the value for the maximum excursions, 
the frequency dependent modified cavitation number, c*(t), 
expressed in equation (7.28) was separated into two parts. 

The first part was the non-oscillating contribution defined 
earlier as a*. The second part was the oscillating contri- 
bution having a magnitude defined by equation (7.29) with a 
frequency equal to the forcing frequency of the shaker table. 


a* = a* + Act* cos wt (7.31) 

The prediction of the magnitude of the pressure oscil- 
lations is given by 


AP = 
s 


P f g 


AG 


sin ft (y/Jt) 
ft cos ft 


(7.32) 


where ft = wA/C. Use of this expression gives rise to the 
frequency varying component Ac* : 

Aa* = 

V £ * 2 


sin ft (y/£) 
cos ft 


(7.33) 


At the entrance plane of the discharge line the ratio 
of y/ A equals unity. Therefore equation (7.33) reduces to 


Aa* 


2gS,G tan ft 



(7.34) 
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The magnitude of A a* will be large when the apparent liquid 
velocity is small. Therefore, for the data obtained in the 
experimental investigation this condition exists for the 
condition of P 0 = 30 psi. This can be explained by the fact 
that the velocity necessary to cause cavitation is lowest 
for the smallest reference pressure. A numerical solution 
for equation (7.34) in the frequency range 

0 F < 70 cps, G = 1 

gives the order of magnitude of A a* noted below 

Ac* * 0.2 (7.35) 

Hence all effects of vibration on the cavitation para- 
meter cause it to be confined to the range: 

3 . 5> < a* ( t) < 3.9 

The limits of possible cavitation numbers are inscribed on 
Figures 26 and 27. It should be noted once again that the 
approximation of one value for the magnitude of the pressure 
oscillation due to the vibration is valid over only a limited 
frequency range. Thus the limits on o*(t) are only inscribed 
over the frequency range where this approximation is valid, 
that is , 


0 < F 60 (cps) 



CAVITATION NUMBER CT (t) s CAVITATION NUMBER 0*10 
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F (CPS) 


26, Time Dependent Modified Cavitation Number Com- 
pared to Approximated Limits for the Conditions, 
d « 1/2, L = 12, T = 80°F , PD = 0, S « 1. 



F (CPS) 

Figure 27, Time Dependent Modified Cavitation Number Com- 
pared to Approximated Limits for the Conditions, 
d = 1/2, L = 24, T - 80°F , PD = 0, S = 1, 



CAVITATION NUMBER ( 7 (t) s CAVITATION NUMBER 



28. Time Dependent Modified Cavitation Number at 

Extended Vibrational Frequencies for Conditions, 
d = 3/8, L = 24, T = 80°F, PD = 180, G = 1, 



Figure 29. Time Dependent Modified Cavitation Number at 

Extended Vibrational Frequencies and High Level 
for Conditions, d = 3/8, L = 24, T » 80°F, PD = 
180, G ~ 5 . 




CAVITATION NUMBER CT (t) £ CAVITATION NUMBER CT (t) 
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F (CPS) 


30. Time Dependent Modified Cavitation Number at Low 
Vibrational Frequencies and High G Level for 
Conditions, d = 3/8, L = 24-, T « 80°F, PD = 0, 

G - 5. 



F (CPS) 

Figure 31. Time Dependent Modified Cavitation Number at Low 
Vibrational Frequencies for Conditions, d - 3/8, 
L = 24, T * 80°F , PD « 0, G * 1. 




CAVITATION NUMBER CT (t) CAVITATION NUMBER <T (t) 
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32. Time Dependent Modified Cavitation Number of Lottf 
Vibrational Frequencies and High G Level for Con- 
ditions, d = 3/8, L « 24, T = 80°F, PD = 180, 
( 1 = 5 . 



Figure 33, 


Time Dependent Modified Cavitation Number at Low 
Vibrational Frequencies and High G Level for Con- 
ditions, d = 3/8, L = 24, T = 80°F, PD = 90, 

G = 5 . 
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7.5 Results and Conclusions for the Cavitation 
Analysis of the Discharge Line 

The use of the modified cavitation number a* defined in 
equation (7.13) shows good agreement with measured data. 

The use of the free streamline analysis in arriving at this 
relation is a useful extension of this theory to flow in the 
entrance of a sharp edged discharge line. Extension of 
equation (7.13) to any similar geometry can evidently be made 
as long as the correct contraction coefficient is known. 

The solution to the cavitation problem in a vibrating 
system was developed considering the liquid to be a hypo- 
thetical fluid, i.e. it can not sustain tensile forces. The 
theoretical significance of the phase deviation, PD, was 
discussed in conjunction with Figure 25. Obviously the ex- 
perimental data presented does not substantiate the importance 
of the phase deviation. 

In order to discuss the apparent lack of sensitivity of 
the liquid in following the pressure oscillations, a brief 
discussion of some of the properties of cavitation that may 
cause this retardation will be made. As mentioned in the 
Literature Survey, Holl (17) showed that the experimentally 
observed phenomenon of the existence of incipient and desinent 
cavitation numbers implies that the process of cavitation is 
more readily induced when the conditions required by the 
desinent cavitation number are satisfied. That is, once the 
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cavitation exists in a liquid, it is easier to maintain the 
condition than to induce it. This fact is in violation of 
the hypothetical fluid concept illustrated in Figure 25. 

The hysterisis effect and a measurable cavitation delay-time 
discussed in the paper by Roll and Treaster (18) further 
implies that the hypothetical fluid model is inaccurate in 
predicting the time dependent response of the fluid. Cavi- 
tation investigations have shown that the onset of cavita- 
tion in a given liquid requires a lower cavitation number 
(i.e. incipient cavitation number, o^) which is manifested 
in a higher liquid velocity than is predicted. Once cavita- 
tion is induced in a system, then the velocity of the liquid 
can be reduced to a lower value while cavitation is main- 
tained. The cavitation number for this lower velocity case 
is called the desinent cavitation number, o^. It does not 
appear that a theory for this cavitation behavior has been 
extablished. 

The work cited in reference (18) yields a possible ex- 
planation other than defining the two cavitation numbers. 

This is the mentioning of a cavitation-delay time for which 
certain general predictions as to velocity dependency, 
dissolved gas content, and surface tension are all parameters 
affecting the onset of cavitation. A significant conclusion 
stated in this reference is that the phenomenon of delay time 
was found to be completely random in nature. Until the 
physical processes are more completely understood the fact 
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that a real liquid may exhibit a delayed response to a loi^ 
pressure ’ must be accepted. Thus the lack of experimental 
phase deviation sensitivity is accounted for by the considera- 
tion that a finite response time exists for the liquid to 
produce sufficient vapor cavities to be visually observable 
as cavitation. An implication is that the pressure fluctua- 
tions associated with high frequency vibrations occur over 
such small intervals of time that the liquid cannot respond 
with a phase change when the static pressure drops below 
the vapor pressure during some part of its oscillation. Of 
course if the mean value of the static pressure drops to the 
vapor pressure, then cavitation could be expected to occur 
continuously. The desinent cavitation number was attributed 
to the fact that cavitation would prevail upon a slight 
reduction in flow velocity which can be considered equivalent 
to a slight increase in the local static pressure. The 
pressure fluctuations attributed to the vibration of the test 
section would be minor enough to not affect the cavitation, 
once started. 

Quantitatively, the limits set on the value of the fluc- 
tuation of A a* seem the best method of interpreting the data. 
Unitl a more detailed analysis is made available, the use of 
the analysis by Schoenhals (44) seems to be a reasonably 
valid approach for a solution to the problem of flow in a 
vibrating discharge line. 
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The fact that the cavitation inception begins at the 
sharp edged entrance to the discharge line has been stated 
with no proof given. Because of the axisymme'tric geometry 
of the discharge line, photographic attempts to record the 
phenomenon were unsuccessful. However \v'ork performed by the 
author in a two-dimensional apparatus allowed a successful 
photographic investigation. Figures 34 are illustrations of 
the two-dimensional apparatus with a geometry very similar ’ 
to that of the sharp edged entrance to the discharge line. 

A cavitation streamer is shown attached to the sharp edge of 
the test section. The same type of phenomenon was visually 
observed in the present experimental work. 
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FIGURE 34a. Stroboscopic Picture of Cavitation 
Along the Free Streamline of a 
Two-Dimensional Discharge 
Line (Magnification 8X) 
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FIGURE 34b. Schematic Diagram Illustrating 
Two-Dimensional System Used in Obtain- 
ing Photograph of Figure 
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CHAPTER 8 

TWO PHASE FLOW IN AN ADIABATIC DISCHARGE LINE AND CHOKING 

8,1 Experimental Observations 
The cavitation discussed in the preceeding chapter is 
the first indication that flashing is about to occur. The 
difference between flashing and cavitation is very subtle 
unless one has had the opportunity to observe both and come 
to a qualitative judgement on differentiation between the 
two. A quote from Collier and Wallis in reference (47) is 
used as an illustration. 

"If saturated liquid at pressure pj flows into 
a region at a lower pressure, p 2 , and is not 
cooled on the way, some of the liquid will 
tend to "flash" to vapor." 

Cavitation has previously been defined as the formation 
of vapor in a liquid when the static pressure reaches the 
vapor pressure because of flow conditions. Cavitation is 
vapor formation due to a local condition of very limited 
extent; whereas, flashing is vapor formation throughout the 
bulk of the fluid. Both phenomena occur due to a pressure 
reduction within the liquid. With this description in mind, 
it is now possible to discuss the two-phase flow conditions 
observed in the discharge line. 
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Upon the initial occurrence of cavitation, several dif- 
ferent regimes of flow existed in the discharge line as 
illustrated in Figure 35 and are described below. 

Regime (1) : Cavitation in the form of vapor 

streamers is observed flowing along the free 
streamline originating at the sharp edged 
discharge line. The vapor streamers consist- 
ing of large numbers of small bubbles dis- • 
appear within the liquid flow once the 
static pres sure' is high enough (con- 
densation) . 

Regime (2) : A vapor annulus forms at the 

entrance and extends in the positive z- 
direction (downward) in response to the 
pressure drop across the discharge line. 

Regime (3): A liquid jet forms from the 

sharp edged entrance with a vapor layer 
between the liquid jet and the wall of 
the discharge line. The jet extends ap- 
proximately tt diameters in the z-direction 
before becoming unstable with the flow 
pattern subsequently transforming to a 
homogeneous one. A normal shock exists at 
the position where the homogeneous flow 
abruptly changes to a single-phase liquid 
flow. 



105 





106 


Regime (4) : This is an extension of the 

previous flow pattern with the exception 
that the shock is now located in the exit 
plane of the discharge line. It cannot 
be driven out of the exit plane in a down- 
stream direction. 

It should be noted that the location of the shock that 
first appears in regime 3 is a function of the exit pressure 
of the test section. Once the shock is established in the 
line, a reduction of the exit pressure forces the shock 
front to move toward the exit plane of the line. The differ- 
ence between flow regimes 3 and 4 described above consists 
O' the fact that in regime 3 a reduction of the exit pressure 
will dislocate the shock while in regime 4 a further reduction 
in the exit pressure does not change the location of the 
shock. The dislocation of the shock front from the exit plane 
of the discharge line into the downstream plenum of the test 
section is physically impossible -when the flow velocity of 
the fluid is reduced to such a low value that the correspond- 
ing static pressure does not allow the existence of a two- 
phase flow mixture. The existence of two-phase flow in the 
discharge line, however, raised the question of maximum flow 
rates, which are achieved under choked flow conditions, for 
given stagnation pressures. 

An investigation of the pressure difference across the 
discharge line as a function of the upstream stagnation 
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pressure, P , showed conclusively that the flow under certain 
conditions could be choked. That is, a maximum flow rate for 
a given constant stagnation pressure value existed beyond 
which any further reduction in downstream pressure did not 
result in an increase of the flow rate. Experimental data 
illustrating this effect are shown in Figure 36 with "best 
fit" curves included. The significance of the figure lies 
in the fact that for a given stagnation pressure, P s , the 
flow rate, Q, reaches a maximum value upon which the pressure 
difference can be increased without a further increase in 
flow rate. The increase in pressure difference at a constant 
upstream stagnation pressure value corresponds to a decrease 
in the exit pressure of the test section. 

8.2 Derivation of an Analytical Model 

The objectives of the study of the problem of two-phase 
flow in an adiabatic discharge line are listed below: 

1. Determination of the maximum flow rate for given 
upstream stagnation conditions and line geometry. 

2. Definition of the pressure losses which occur in 
this type of flow. 

3. Determination of the parameters that directly affect 
the normal shock which appears in the line. 

The first objective noted above for a two-phase flow 
situation can best be achieved if a comparison with choking of 
a perfect gas in a converging-diverging nozzle is made. 
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Figure 36. Volumetric Flow Rate 'f *o V Various Upstream 

Stagnation Pressure Conditions as aFun'ction of 
Pressure Drop Across the Discharge Line. 
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Consider a converging-diverging nozzle with constant 
upstream stagnation conditions and a variable exit pressure. 
The working fluid is an ideal gas. When the exit pressure 
equals the supply pressure, the trivial case of no flow 
prevails. As the exit pressure is lowered, the nozzle 
behaves like a typical venturi tube. The flow rate increases 
and the static pressure in the throat becomes sensitive to 
the exit pressure. This’ type of flow behavior continues 
until an exit pressure is reached where the flow rate at the 
throat becomes sonic. From this value of the exit pressure 
on to all lower values the flow rate stays constant. This 
operating condition is defined as choked flow. A further 
increase in flow rate would require an increase in upstream 
stagnation pressure or a change in temperature. The only 
effect that a change in the exit pressure produces is a 
change in the location of the shock which forms in the diverg 
ing portion of the nozzle. 

The f 1 ow of a liquid in a discharge line behaves in a 
manner very similar to that described above for a perfect gas 
That is, the liquid flow rate at a given upstream stagnation 
temperature and pressure depends on the exit pressure in a 
given system. Once a minimum value of the exit pressure is 
reached no further increase in flow rate can occur without 
changing the upstream conditions. This minimum pressure that 
the liquid can sustain at a given fixed temperature is its 
vapor pressure. Consideration of a'Control Volume having 
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its exit plane in the two-phase mixture shows that the maxi- 
mum flow rate will occur once the exit pressure is identical 
to the vapor pressure of the liquid. An attempt to further 
reduce the exit pressure of the liquid would only result in 
phase change of the liquid and hence an increase in the 
volumetric flow rate; however, the mass flow rate would 
remain constant. Experimentally this type of investigation 
would remain constant. Experimentally this type of investi- 
gation would be extremely difficult to perform because it 
would either require a pump which can handle a two -phase 
mixture or a very large tank maintained at constant pressure. 

Using this criterion for the discharge line problem, 
the application of the free streamline solution presented by 
Rouse and Abul-Fetouh (6) provides the limiting flow condi- 
tions. Once the pressure along the free streamline reaches 
the vapor pressure and the pressure drop along the discharge 
line is sufficient for the jet to develop fully, the static 
pressure in the cylindrical cross section of the discharge 
line will be constant with a value equal to the vapor pressure. 
Since no further reduction in the pressure is possible, the 
condition for maximum flow now exists. It has to be noted 
that when the pressure in the liquid jet reaches the vapor 
pressure, nucleation can occur and the character of the flow 
changes. The flow in the discharge line has been observed to 
change from the jet flow to a homogeneous flow pattern with 
a normal shock present. 
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The defining relation for choked flow in a discharge 
line with a sharp edged entrance at a fixed stagnation 
pressure and temperature is given below. It is a limiting 
condition for the application of Bernoulli's equation along 
a streamline. 


or 


P v + 


1 

2 Pf 



( 8 . 1 ) 



( 8 . 2 ) 


Data for several geometric configurations and for different 
temperatures are illustrated together with the solution to 
equation (8.2) in Figure 37. 

The theoretical curves in Figure 37 represent the maxi- 
mum possible values of the apparent liquid velocity corres- 
ponding to given values of the stagnation pressure, i.e, 
they are by definition the choking curves. For a fixed value 
of the upstream stagnation pressure, , all values of the 
apparent liquid velocity up to and including the value rep- 
resented by the solution to equation (8.2) are possible. 

In Figure 37, these possible flow conditions correspond to 
the areas above the curves. Since these choking curves are 
plotted as functions of the apparent liquid velocity, a 
straightforward* conversion of V^* to the mass flux is possi- 
ble by simply multiplying values of V f * by the density of 
the liquid. In this manner the mass flow rate is readily 
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Figure 37. Measured Flow Rates Under Choked Conditions as a 

Function of Stagnation Pressure for Various Discharge 
Line Geometries Compared with the Analytically Pre- 
dicted Maximum Flow Rates. 
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determined by knowing the upstream stagnation conditions and 
the geometry of the discharge line under investigation. 

Since the maximum flow rates (choked flow) are deter- 
mined by equation (8.2), the analysis must proceed in a more 
specific, manner to investigate the flow within the discharge 
line in view of prevailing static pressures and temperatures. 
Of primary importance is the prediction of the static pres- 
sures and pressure losses incurred by the fluid flowing 
through the discharge line. An explanation of the notation 
and geometry used in the subsequent derivation is illustrated 
in Figure 38. 

The static pressure measured along the wall of the dis- 
charge line in the region 

0 < z < p 

equals the vapor pressure of the liquid, and is identical to 
the static pressure along the free streamline. The liquid 
jet, however, maintains a pressure field which is a function 
of the coordinates r, 0, z until its width is close to its 
asymptotic value. Within the jet the static pressure is not 
constant until z ~ x\. 

For 0 < z < q 

^static,jet ~ z ) (8.3) 

This functional relation is graphically illustrated by Rouse 
and Abul-Fetouh (6), At the position where the static pressure 
is constant in a radial direction, the jet changes into a 
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LIQUID FLOW 



'H = LENGTH AT JET TRANSFORMATION 

(HOMOGENEOUS FLOW STARTING POlfoT) 

f = POSITION OF TWO-PHASE * SINGLE PHASE 
SHOCK 


L * LENGTH OF DISCHARGE LINE 

Figure 38. Geometry and Coordinate System for Analysis of 
Adiabatic Two-Phase Flow in Discharge Line. 
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homogeneous, two-phase flow pattern. In the homogeneous 
flow regime, the pressure is constant and equal to the 
vapor pressure. Experimental verification of this fact is 
presented in Figure 39 where the pressure difference between 
the upstream stagnation pressure and the static pressure, 
p^, is illustrated. The location of the pressure tap through 
which the static pressures, p 2 , were measured is shown in 
Figure 7. The choked flow curve shown in Figure 39 is 
derived from the following expression. 

P s ‘ P 2 = 7 Pf (V £ */C c ) 2 (8.4) 

with p 2 = P y 

For comparison purposes, a similar set of experiments 
was performed where the flow conditions were adjusted in such 
a way that no cavitation or flashing occurred. Again the 
static pressures in the discharge line were measured and for 
this situation also displayed in Figure 39. The difference 
between the given stagnation pressure and the measured static 
pressure, p 2 , is shown as the ordinate while the apparent 
liquid velocity, V^* , which, now in the case of liquid flow 
only, is the real one-dimensional velocity, appears as the 
abscissa. 

If in a given experiment with constant stagnation pres- 
sure, say P s = 40 psia, the exit pressure is decreased, i.e. 
the flow velocity, V^* , is increased', then the difference 
between the constant stagnation pressure and the measured 
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Figure 39. Comparison of the Dynamic Pressure for Single 
Phase Liquid Flow as a Function of Apparent 
Liquid Velocity with Difference of Stagnation 
Pressure and Static Pressure of a Flashing Fluid 
for Various Stagnation Pressures. 
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static pressure follows the lower curve as predicted by equa- 
tion 6.5 shown in Figure 39. Once flashing occurs the static 
pressure, p^ , reaches the vapor pressure of the liquid and 
the difference between the stagnation pressure and the new 
static pressure suddenly increases to a higher value. These 
observations are shown for various constant stagnation pres- 
sures [dashed lines) in Figure 39. Since the pressure tap 2 
was located a sufficient distance from the entrance to the 
discharge line, the flow was always choked when the measured 

static pressure was equal to the vapor pressure. 

\ 

The static pressure in the two-phase flow region of the 
discharge line was found to be equal to the vapor pressure. 

Now consideration of the stagnation pressure losses must be 
made. The pressure loss expression for the sudden contraction 
at the entrance defined for the single-phase liquid in equa- 
tion (5.2) can be modified and applied to the contraction 
with flashing. Richardson (49) found a value for the loss 
coefficient K to be approximately 0.20 for the area ratios 
used in the present test section configuration. It should 
be noted that Richardson dealt with two-phase, two-component 
flow; he did not study flashing. Richardson's value for the 
loss coefficient was defined on the basis that the fluid 
density equaled the mixture density rather than that of the 
single-phase liquid. This implies that the void fraction 
must be included in the definition of the entrance loss. 

That is , 
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where K = 0.20 

As demonstrated later, this assumption for the loss at the 
entrance is justified by the experimental results obtained 
in this investigation. 

The entrance loss defined in equation (8.5) accounts 
for the stagnation pressure loss due to the contraction and 
jet formation. This definition explains the pressure loss 
from the entrance plane of the discharge line to the position 
n. In this length of the discharge line, the liquid phase 
of the two-phase mixture does not come into contact with the 
wall. The vapor annulus separates the liquid jet from contact 
with the wall, and subsequently the frictional losses over 
this portion of the flow are negligibly small. This condi- 
tion is stated below in a more explicit form. 


For 0 < z < n 

iP friction ,TP 


( 8 . 6 ) 


However, when the transition to homogeneous flow takes 
place in the discharge line, a significant pressure loss due 
to frictional effects does occur. The friction factors ob- 
tained previously for the single-phase flow were applied by 
taking the two-phase friction factor, f Tp , equal to the already 
known single-phase liquid friction factor. This is illustrated 
by Collier and Wallis (51) who discuss a suitably defined 
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two-phase flow friction factor for the homogeneous model, 
i.e. pseudo-fluid. 

f TP « f * * E/d] (8.7) 

f 

where e/d = relative roughness 

The frictional pressure loss can be defined as follows. 

For n < z < £ 


l£ 
3 z 


f TP 

“ 3 ” 


ip V : 
2 p m m 


( 8 . 8 ) 


i friction ,TP 

Obviously the void fraction, a, needs to be known before 
equation (8.8) can be used because of the functional depend- 
ency of the density and velocity of the mixture on a. 

Consider now the flow at two locations in the discharge 


line defined below. 

z * - n - dz 

1 . 

z ! = n (8.9) 

2 

The continuity and momentum equations are written and simpli- 
fied between the two control planes. 


Continuity: G x ' = G 2 ' (8.10) 

Momentum: 

P lt + (G f V £ + G g V g ) t , « p 2 < + (G f V £ + G g V g ) a , 

( 8 . i i ) 

Note that over- the interval from z l , to z 2 , the frictional 
force in the z-direction is still defined as being zero. 

Here the liquid phase does not come into contact with the 
wall until z = n; thus the frictional pressure gradient can 
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be neglected. The static pressure in the flashing flow 
regime has been shown previously to be equal to the vapor 
pressure . 

That is, 

Pi ' = P 2 ' = p v (8.12) 

Since the vapor in the region 0 < z < tj flashes from the 
surface of the free streamline envelope, the component of 
velocity of the vapor in the positive z-dircction can be 
assumed to be equal to the free streamline velocity. 


V fi ' = V • (8.13) 

With the assumption of homogeneous flow at , , another con- 
dition can be stated. 


V f 2 ' = V g 2 ' = V m (8 .14) 

Substitution of equations (8.12 - 8.14) into (8.11) yields a 
modified momentum equation. 

( G f * G g),. V fl , - (C f + G g ) 2 , V m (8. IS) 

From continuity, 

G i’ = ( G f + G g )i» = G 2 . = (G f + G g j 2 , (8.16) 

it is obvious that 

The theory of jet flow claims that the area of the jet equals 
the area of the orifice times the contraction coefficient. 

The velocity of the liquid at n is defined below. 
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Qr 


•r 'x V £ 

V £,, = A]7t ’ C c CndVj|) = C7 


(8.18) 


With the assumption that the density of the vapor phase is 
much smaller than that of the liquid phase, 

Pg << P f (8.19) 

the continuity relation given in equation (8.10) is rewritten 
as 

''m 'm 


p m V m = p f V f * 


( 8 . 20 ) 


since 


P m = ap g + (1 - a) p f = (1 - a) p f 


V « 


v f* 


( 8 . 21 ) 


m 1 - a 

Upon combining equations (8.17), (8,18), and (8.21), an ex- 
pression for the value of the void fraction at z => n is 
obtained , 


V 


v * 

v f 

1 - a 


( 8 . 22 ) 


Therefore 


a s 1 * C c at z s n (8.23) 

With a known value for the void fraction, a, in the 
plane of transition from jet flow (reversed annular) to 
homogeneous two-phase flow, it is now possible to analyze 
the shock phenomena observed in the discharge line under 
choked flow conditions. 
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Consider the situation shown in Figure 38 where the 
normal shock is located at a position £ in the discharge 
line. One-dimensional normal shock relations describing 
the flow dynamics can be used to derive the downstream flow 
conditions in the single-phase liquid. 


According to Figure 38 the coordinates are written 


below. 

I,et z = £ - dz/2 

( 8 . 2 * 

z 2 = t + dz/2 

At location z l , two-phase homogeneous flow exists. At loca- 
tion z 2 , single-phase liquid flow prevails. Treating the 
homogeneous regime as flow of a pseudo-fluid allows the fol- 
lowing relations to be formulated. 


J = V 
2 . 2 


J = V 
1 1 


P, ~ Pf 


p i p m 


Vr = V 
f 1 gl 


where' J ~ volumetric flux at position 2, etc. 

The continuity equation is applied for the control areas at 


z, and z 
1 2 


P X \ " P , " p 2 V 2 = P 2 J 2 " P f V 


(8.26) 


or J = (8.27) 

1 p ! 

The volumetric flux difference (relative flux) is determined 
to be 

J 2 - J i = J J 1 " p f/°x) (8.28) 
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Pf v* 

J a - (1 - p £ / Pl ) 


(8.29) 


With the assumption that 


P x = C 1 " a) p£ 


it becomes 


J - J = 

2 l 


Pf v 


1 - 


1 - a 


(8.30) 


which simplifies to 


J - J = 


-a V £ * 


2 i 1 - a 


(8.31) 


The momentum equation written for the coordinates z y and z 2 
can be expressed in the form below. 


p + GJ = p + GJ (8.32) 

1 12 2 

Once again the fact that the pressure in the flashing flow 
regime equals the vapor pressure is used to simplify equation 
(8.32) . 

P 2 - P y - G (J 2 - J x ) (8.33) 

Substitution of equation (8.31), the continuity, equation, into 
equation (8.33), the momentum equation, yields the expression 
for the pressure in the liquid immediately downstream of the 
normal shock. 


P 2 


P v + I p f v f‘ : 


2.a] 

1 - a 

\ 


C 


(8.34) 
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From the position z 2 to the exit plane of the discharge 
line, the one dimensional friction loss relation for a 
single phase liquid, equation (6.3), is applied. Again the 
expansion loss defined in equation (6.6) also is applied at 
the exit plane since the flow is single phase after passing 
through the shock. 

Use of the energy equation between coordinates z 2 and 
the exit region of the test section determines the downstream 
stagnation pressure in the test section, P e . 


P e * P 2 + P£ g (L - o * | p f V f * 2 


AP 


friction 


L 

- AP 

^ expansion 


( 8 . 35 ) 

Substituting the relations for the friction and expansion 
losses in addition to the definition of pressure p 2 into equa- 
tion (8.35) yields the somewhat complicated expression. 


p +i„ V * 2 
v 2 p f v f 


2a 


1 - a 


g (L - C) 


‘£ V f 


ilk - €) 1. o 
a i p 


v * 

f v f 


- (1 - d 2 /D 2 ) 2 \ p f V f * 2 


(8.36) 


After rearranging terms and simplifying, the exit stagnation 
pressure is given by equation (8.37). 


P 

e 


1 

J 


p £ 


V 






+ \ + P f g (L - 


(8.37) 
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The total stagnation pressure loss across the discharge line 
due to fluid flow can now be determined. Note that the losses 
are zero for the condition of no flow in this definition. 

AP = P s - P e + P f g L (8.38) 

If flow is occurring within the discharge line, then the 
pressure losses due to the flow must be subtracted from the 
exit stagnation pressure, P . Thus for flow in the discharge 
1 ine , 

AP = P s + p f g L 

e ’ ^entrance $ n cTz f jp + S P m ( z )d z 

’ (8.39) 

Substitution of equation (8.37) into equation (8.39) yields 

the expanded relation for the stagnation pressure loss caused 

by flow in the discharge line. 

AP = P s - P y + p f g L + Ap entrance 

♦ I p£ ?) + u - 

dz - T g 0 m Cz)d2 - p f g(L - E) (8.40) 

To simplify equation (8.40), it is possible to combine the 
hydrostatic pressure terms into one simplified relation. 

£ £ 

p f g L ' J g P m ^ z ^ dz " Pf g ( L * £) = / 0 g (Pf " P m ( z )J dz 

(8.41) 

Substituting equations (8.5) and (8.41) into the pressure loss 
equation (8.40) with the assumption, 


f 9P l 

f ,TP 
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for 0 ± z £. E, 


[P f ' p m ( z )l = to(z) p f l' 


(8.42) 


yields the final form of the pressure loss equation. 

A p = \ P f V f * 2 |o.2(l- a ) j + UlfZl + (l-d 2 /D 2 ) 2 - i~j 1 

n r a 


+ p 


p v * r„ [H|f jXP dz M; s Pf *C*)dz (8.43) 


With equation (8.43), it is now possible to predict the 
loss in - recoverable stagnation pressure across the test 
section for a given position of the shock. This loss is 
measured by the differential pressure transducer located 
between the inlet pipe and exit pipe of the test section. 

It was adjusted to read zero under the no flow condition, 

V^* =0. It will be demonstrated later that equation (8.43) 
can be used for the prediction of the location of the normal 
shock front under given flow conditions. 


8.3 Discussion of Choking Conditions 
In this section, the definition for choking in two- 
phase flow is applied to the analytical model described in 
the previous section. This is done to illustrate that the 
defining equations satisfy the accepted choking conditions. 
Use is made of equation (8.43) for the analytical prediction 
of choking. These predictions are then compared with experi- 
mental measurements made on the discharge line. 
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A one -dimensional steady state analysis of homogeneous 
two phase flow using the momentum equation yields an expres- 
sion for the pressure gradient in a system which is the sum 
of several components, each of which is associated with a 
particular physical effect. The homogeneous flow model is 
employed since this flow regime is clearly' discernable in 
the discharge line upstream of the shock, and because it is 
analytically the most tractable to apply. In analytical form 
the pressure gradient in the positive z-direction is given by 


where 



(8.44) 



frictional component of pressure gradient 



acceleration component of pressure gradient 



gravitational component of pressure gradient 


Equation (8.44) can be found in the literature in a 
multitude of forms in terms of various .combinations of vari- 
ables. A pertinant discussion is presented in references 
(S3) and (54). To illustrate the significance of equation 
(8.44), it can be written in a steady-state analysis as a 
combination of influence coefficients. 


C £ ' 


friction ■ 

c x 

Influence coef- 

phase change 

I 

* 

ficient of 

area change 

c g 1 

i 


gravity 

t * ' 


on pressure 
gradient 
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Using the concept of an influence coefficient, equation 
(8.44) can be shown to have the following form. 



P i r U-X 1 y-1 1 dA . r-1 

c f C x 3Tz C A X Jz C g g cos T 


(8.45) 


The term, M z , in the denominator represents the Mach number 
squared of a pseudo fluid, and it has the same significance 
as in the compressible single component flow equation. If 
an. analogy is made between single phase compressible fluid 
£low r and the homogeneous "pseudo" fluid, a choking criterion 
can be established. For compressible flow, a condition of 
Mach number equal to unity at some location along the flow 
path is always associated with choking. The' above derivative, 
equation (8.45), becomes infinite in value. For homogeneous 
flow, this derivative is indeterminate not only. for Mach 
equal to one but also for sudden changes in the quality, X. 

Consider now' the flow at position £ indicated in Figure 
38. At this location a finite change in quality occurs in 
zero length at the interface between the homogeneous "pseudo" 
fluid and the single phase liquid. The relation, equation 
(8.45), becomes infinite at this location; and although the 
Mach number does not necessarily equal unity the flow is choked. 

A special case of the flow in the discharge line can now 
be discussed. Let the shock be located immediately downstream • 
of the transition from jet flow to homogeneous flow. The 
following conditions therefore exist. 
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(e - n) * 0 


(8.46) 



(8.47) 


Following from the above assumption, the relation for the 


void fraction can be determined from the contraction coeffi- 


cient . 

. •• . 2 - C 

1 + ot [ _ 1 + a i _ c 

1 - a 1 ^ " 1 - a' n C c 


(8.48) 


Substituting the above conditions into the expression 
for the pressure loss across the discharge line, equation 
(8.43), yields a relation which represents the minimum pres- 
sure loss for a choked condition. Neglecting the contribu- 
tion of the hydrostatic pressure change results in 


J g ct(z) p f dz ~ 0 


(8.49) 


Therefore 


■W - 7 Pf V 2 [ 0 - 2 u-«V £i r il * ' 1 - dJ/D2)2 ' t; 


+ P s ' P v CS.50) 

Once again an analogy can be made between the flow in 
the discharge line and single phase compressible flow. In 
particular, the pressure profiles for a single phase compres- 
sible flow in a converging-diverging nozzle serve as a basis 
for comparison. There exists a back pressure for the nozzle 
at which the flow reaches a sonic condition in the throat 
which can be considered as an infinitely weak shock, and then 
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exits subsonically with the diverging portion o£ the nozzle 

l 

acting as a diffuser. At a value of the back pressure loiter 
than this "critical*' pressure a normal shock exists in the 
nozzle, the location of which is dependent on the value of 
the back pressure. A similar situation prevails in the dis- 
charge line. There exists a back pressure where an infinite- 
ly weak shock exists in the liquid. This occurs under the 
condition of ri = £. A further decrease in the back pressure 
below this "critical" pressure moves the shock in a down- 
stream direction without affecting the flow rate. 

Equation (8.50) can be used *o determine the maximum 
value for the pressure loss for non-choked flow. A compari- 
son with experimentally measured values for the pressure loss 
in two geometric configurations where the shock was located 
in what is usually the transition region from jet flow to 
homogeneous flow, £ - n ~ 0, is given in Table 4. 

Table 4 


Comparison of Experimental and Theoretical Values 
of the Minimum Pressure Loss for Choked Flow 

(€ - n) * 0 



L = 24, 

d = 0.375 

, f = 0.0155, T 

= 80° F 

p s 


V f * 

... AP .t.h.e ore t i c a 1. „ 

^experimental 

68 


61.2 

61.2 

61 ' “ 

78 

; 

65.6 

70.0 

68 

87 

* 

69 . 3 

78.5 

77 


L = 12, 

d = 0.5, 

f = 0.0160, T = 

H«i 

O 

o 

oo 

p s 


V 

^theoretical 

A p exper imental 

28 


39.0 

18.6 

18.5 

32 


41.8 

21.0 

22 

42 


48.0 

28.0 

29 

50 


52.4 

33.5 

34 
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The special case shown in Table 4 has already been 
described as the minimum pressure loss across the .discharge 
line for choked flow. This means that “the exit pressure 
could continue to decrease for a constant upstream stagna- 
tion pressure without increasing the flow rate. The only 
effect of such a reduction would be to move the shock towards 
the exit plane of the discharge line. The functional rela- 
tionship which exists between the motion of the shock and the 
exit pressure can be determined by a differentiation of equa- 
tion (8.37) which relates the exit pressure to the shock 
coordinate, £ • 

P ff) = P 
e'- w e 


£=n 


3 AP 

as 


(S - n) 


(8.51) 


The expression for the stagnation pressure loss, equa- 
tion (8.43), can be differentiated with respect to the co- 
ordinate That is, 


3(6P) . L, v 

2 P f\f 


3?{0-2U-<»)l n * * Cl-dVD 2 ) 2 - i^| } 


* TZ f l}t}f,TP Jz + h s p f a ( z > dz 


S' 

(8.52) 


Simplifying this relation gives, 

3 ( AP) f 1 „ „ * 2 ^ a f3p' 

TT 1 S' 7 p f v f + 3T \ ISz; 


f ,TP 


dz 


3T 


L g Pr a(z)dz 


(8.53) 


The definition of the frictional pressure gradient was 
given in equation (8.8) in terms of the mixture density and 
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mixture velocity. The continuity equation can be used to 
redefine the frictional pressure gradient in terms of the 
void fraction and the apparent liquid velocity. This form 
of the gradient is shown below as equation (8.54). 


3p 


9 Z 


TP 

“cT 


2 p f V 


f,TP a U " «) 

The following equation is obtained by substituting 
equation (8.54) into equation (8.53). 


(8.54) 


3 UP) 


f 1 


V, 




.2 , ■ l TP 1 

~T~ 2 p f 


v * 2 L_ f‘ 
V f 35 J 


1 - a(z) 


dz 


3C 


/ g P f <x(z) dz 


(8.55) 


Combination of terms in equation (8.55) results in a simpli- 
fied expression for the change of pressure loss with respect 
to a change in the coordinate 5* 


3 UP) 


f 1 

= a* 2 


Of Vf*' 



1 


+ ||- g p f a(z) dz (8.56) 

To use equation (8.56) the functional relationship of 
a with respect to z must be specified. Under the assumption 
that the void fraction undergoes only very minor changes in 
the flashing mixture as a function of displacement along the 
discharge line, the following statement can be made. 
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for n < z < 5 

a(z) - a(n) = constant = 1 - C (S . 57) 

o 

Applying (8,57) to the exact relation, equation (8.56), 
a much simplified expression for the pressure loss dependency 
is obtained,' 

+ g P f a (8.58) 

h 

If the equation (8.58) is applied to high speed flows, the 
apparent dynamic pressure term is very much larger than the 
hydrostatic pressure term. Thus the additional simplifica- 
tion given as equation (8.59) can be applied to (8.58). 

That is, with 

* « 7 Pf V* J [rM I (8.59) 

the final form of ths pressure dependency upon the coordinate 
is obtained, 

d( AP) --1 it ** i f a i! 

3? I p f V f' a - (r ” a] (8.60) 

A comparison of the numerical results obtained from 
equation (8.60) with experimental data is shown in Figures 
40 and 41. Since this equation represents the derivative 
(i.e. slope) of the pressure loss with respect to 5 , it is 
illustrated in the Figures 40 and 41 gives lines having slopes 
equal to the values obtained from the solutions of equation 
(8.60). The graphical representations of equation (8.60) are 
inscribed in these figures in positions which most closely 
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40 


EQUATION (8.60) 

TABLE 4 VALUES ® 
(MINIMUM PRESSURE 
LOS: 




d= 0.5 


£ (INCHES) 

Figure 40, Stagnation Pressure Loss in a Discharge Line of 

d - 0,5 and L ~ 12 as a Function of Shock 'Position 
'"J'ith Various Upstream Stagnation Pressures, 



AP (PSD 



6 (INCHES) 

Figure 41. Stagnation Pressure Loss <in a Discharge Line of 
d = 0.375 and L = 24 as a Function of Shock Posi 
tion with Various Upstream Stagnation Pressures. 
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Figure 42. Minimum Pressure Losses as Functions of Apparent 
Liquid Velocities, i.e. Upstream Stagnation Pres- 
sures, for Various Geometric Configurations of 
the Discharge Line. 



13 8 





V f (Ft/sec.) 



Figure 42 
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match the experimental data. Also included in the figures 
are the theoretical minimum pressure losses noted in 
Table 4, 

A final comparison of the analytical pressure loss 
relation is made by consideration of the minimum pressure 
loss for choked flow as expressed by equation (8.50). Since 
it was defined as a minimum loss relation, all experimental 
pressure loss data for choked flot^ necessarily must have 
values greater than or equal to that predicted by the solu- 
tions of this equation. Data for several different geometric 
configurations of the discharge line are compared with the 
theoretical minimum loss value in Figure 42. To prove the 
validity of equation (8.50)', all data must lie on or above 
the line representing the solution to (8.50), which is 
indeed observed. 

8.4 Investigation of Choking Conditions 
for a Vibrating Discharge Line 

In the preceeding sections, the theory of two-phase 
single-component non-vibratory adiabatic discharge flow was 
developed. An extension of this theory to a vibrating dis- 
charge system is discussed in this section. Primary con- 
sideration is given to the effect which vibration has on two 
phase flow choking. 

The fact that the flow rate attains a maximum value 
when the jet in the entrance region transforms into a 
homogeneous two-phase flow with a shock present has been 
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discussed. For the flow in a vibrating discharge line, the 
static pressure has an oscillating component super-imposed 
upon its mean value as shown by Schoenhals (44) and illus- 
trated in this study in equation (7.32). The fact that the 
sinusoidal variation of the static pressure can only be 
sustained in the liquid phase becomes obvious when the 
concept of a hypothetical fluid is used. Once again, the 
vapor can only exist at one pressure value as long as suf- 
ficient liquid is available to undergo a phase transforma- 
tion to the vapor phase and the vapor maintains an interface 
with the liquid. If a pressure pulse were to propagate 
through the single-phase liquid to a liquid-vapor interface, 
the pressure wave would be absorbed by phase transformation 
processes. Whether this pressure wave is a compression wave 
or a rarefaction wave does not matter. The conditions of 
thermodynamic equilibrium require that the two-phase mixture 
will be forced to have sufficient amounts of vapor condense 
to liquid to maintain the mixture pressure at the vapor 
pressure of the liquid for a compression wave, and the two- 
phase mixture will be forced to let the liquid phase evapor- 
ate in sufficient amounts in order to maintain the mixture 
pressure at the t'apor pressure for a rarefaction wave. 

Consider now the case of a real fluid for which transi- 
ent effects have to be included. This is to say that the 
fluid cannot respond immediately to a change of pressure. 

If the same pressure wave were to propagate through the 
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single-phase liquid to the liquid-vapor interface; the 
vapor could possibly sustain the wave for a brief period 
until equilibrium could be reestablished. The real two- 
phase liquid-vapor fluid would initially damp the pressure 
wave and finally absorb it if sufficient fluid were availa- 
ble. The amount of two-phase mixture necessary to absorb 
a wave of known strength would of course depend on the rate 
processes of the particular fluid involved and such consi- 
derations as interfacial area and impurities. 

Applying the concepts stated in this brief discussion 
to the test fluid as it passes into the discharge line, the 
ratio of interfacial area to liquid volume can be used as a 
criterion for the rate of maximum phase transformation. The 
region of greatest interfacial area per unit volume would 
exist in the homogeneous two-phase flow mixture. 

An experimental check on the fluid's behavior when 
choked was performed. The high frequency pressure transducer 
measuring pressure p t in Figure 6 showed that when the vapor 
annulus existed between the liquid jet and the wall of the 
discharge line, no oscillating component of the static pres- 
sure was detected. This implies that the test fluid was 
capable of absorbing pressure oscillations when the two-phase 
liquid-vapor combination existed. 

Because the static pressure in the two-phase region of 
the discharge line under vibrating conditions did not differ 
from that which existed in absence of vibration, no change 
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Figure 43. Measured Flow Rates Under Choked Conditions in 

Presence of Vibration as a Function of Stagnation 
Pressure for Various Discharge Line Geometries 
Compared with the Analytically Predicted Maximum 
Flow Rates in Absence of Vibration. 
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in flowrate for given upstream stagnation conditions was ex- 
pected for the case of longitudinal vibration compared to 
the case where vibration was absent. Data taken for various 
geometries, acceleration levels, and frequencies are comparec 
to the prediction for the non-vibrating discharge line, equa- 
tion (8.2), in Figure 43. No detectable difference from the 
theoretical prediction was noticed. 

8.5 Results and Conclusions for the 
Two Phase Flow Analysis and Choking in the Discharge Line 

Application of the free streamline theory to the flow 
pattern in the entrance of a sharp edged discharge line 
yielded a solution which was used to predict the maximum 
rates under choked flow conditions. Comparison of the theo- 
retical prediction with experimental data in Figure 36 illus- 
trates the accuracy of the prediction. It should be noted 
that while experimental values of the parameter, d 2 /D 2 , re- 
stricted the correlation between data and theory to the 
special case of the limiting value of the contraction coef- 
ficient, C c = 0.611, the analysis can be applied to any 
similar geometric arrangement as long as the proper value 
of the contraction coefficient is available. 

The advantage of this type of analysis was realized in 
transforming the multidimensional flow problem in the 
entrance region to one which could effectively be treated 
as a one dimensional problem. This technique consisted of 





an extension of the single phase liquid flow analysis pre- 
sented earlier to a two-phase flow analysis. 

The entrance loss coefficient for the single phase 
liquid flow was found -to be very much greater than that of 
the two-phase flashing flow. The modified correlation from 
Richardson agreed well with experimental data, although its 
original application was directed towards a two -component, 
two phase flow system. As in the case of the single phase 
liquid flow, the greatest pressure loss occurred in the ex- 
pansion process rather than in the contraction process. 

Since the single component flashing flow in the discharge 
line represented a two phase flow during the transformation 
from jet flow to homogeneous flow, the similarity between 
Richardson's flow problem and this flow situation was suffi- 
cient to allow the application of a similar relation. In 
the entrance loss concept frictional losses were disregarded 
until the jet expansion forced the liquid into contact with 
the wall. 

During the transformation from jet flow to homogeneous 
flow, the mixture velocity was shown to be equal to the jet 
/elocity. At the coordinate n , this velocity equaled the 
/elocity along the free streamline as formulated in the 
analytical solution. The existence of continuous liquid 
velocities in the region of flow transformation allowed a 
prediction for the void fraction at the point of inception 
of homogeneous flow. The void fraction at this point is 
dependent upon the contraction coefficient. 
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The determination of the void fraction at the position 
n made it possible to predict the minimum pressure loss for 
choked flow conditions. Minimum pressure losses occurred 
when an infinitely weak shock existed in the region of trans- 
formation from jet to homogeneous flow. The correspondence 
between the theoretical prediction and the experimental data 
shown in Table 4 for this special case is extremely good. 

An extension of the pressure loss relation to predict 
the location of the normal shock (i.e. coordinate 5), was 
presented. With the assumption of a constant void fraction, 
the expression derived in the analysis gave the shock posi- 
tion dependency as a linear function of the' pressure loss at 
a given flow rate. 

Two photographs of the normal shock taken with the aid 
of a stroboscopic flash unit are presented in Figures 44 and 
45. The duration of the flash was approximately 9p seconds. 
The dark region upstream of the shock interface is the homo- 
geneous mixture. The numerous interfacial regions between 
liquid and vapor provided many areas of total reflection 
which effectively reduced the amount of transmitted light 
energy. The lighter area downstream of the interface is the 
single phase liquid. Depth of field of the lens and the 
cylindrical cross section of the discharge line prohibited 
sharper photographs. However, the contrast of the interface 
is still apparent. 
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Data from the experimental investigation into choking 
of the discharge line under the influence of longitudinal 
vibration were compared with the analytical prediction form- 
ulated for the non-vibrating line. The excellent agreement 
between the data and predictions proved that the vibration 
did not affect the maximum flow rates for given upstream 
conditions of the liquid. No determination of the pressure 
losses was attempted for the vibrational cases. 
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FIGURE 44 . Stroboscopic Photograph of Normal Shock 
in Discharge Line (Magnification 3X ) 
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FIGURE 45 . Stroboscopic Photograph of Normal Shock 
in Discharge Line (Magnification 4X ) 
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CHAPTER '9- 
SUMMAPY 

This study was concerned with both single phase liquid 
and two-phase liquid-vapor fluid flow in a discharge line 
with a sharp edged entrance. Initial efforts were spent 
determining the applicability of the classical single phase 
liquid flow relations to this particular type of problem. 
Following this initial phase of the overall study, a cavita- 
tion investigation was performed with emphasis on determin- 
ing a suitable cavitation parameter which could be used to 
predict the onset of two-phase conditions within the discharge 
line. An extension to totally two-phase flow conditions in 
the discharge line followed the cavitation work. The totally 
two phase flow condition was attributed to flashing of the 
liquid to its vapor phase because of the existing flow condi- 
tions. All of the flow problems were investigated both in 
the absence and in the presence of longitudinal vibration of 
the discharge line. 

In order to obtain experimental data for comparison with 
analytical predictions, a two-phase fluid flow loop facility 
having a test section mounted on an electrodynamic shaker was 
designed and constructed. The entire experimental system was 
a combination of the flow loop and the shaker facility. This 
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system was quite flexible in terms of operating regimes having 
a fluid temperature range from 60° F to 180° F, flow rates 
from 0 to 45 GPM, vibrational frequencies in the range of 
10 to 2000 cps and maximum acceleration levels up to 20 g 
for the test sections used on the shaker. 

The check-out of the ex*per imental facility was done 

/ 

with single phase liquid- flow through the test section under 
non-vibratory conditions. The instrumentation of the facility 
was checked for correct operation by using well known analyti- 
cal relations of single phase liquid flow as a basis for com- 
parison with experimental data. 

More specifically, the one-dimensional analysis presented 
in Chapter 6 describing the stagnation pressure loss across 
the sharp edged discharge line for single phase liquid flow 
yielded an expression the solution of which compared favorably 
to experimentally obtained data. An analytical solution to 
the pressure loss across the discharge line under the influence 
of longitudinal vibration was not attempted; however, an experi- 
mental investigation was conducted. The results of this work, 
presented at the end of Chapter 6 , showed in general that the 
presence of this type of vibration increased the existing 
pressure losses. A trend of larger losses for lower frequen- 
cies at a given acceleration level was noted. Also, at a 
given fixed frequency, larger acceleration levels usually 
produced larger pressure losses. 
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By treating the flow at the sharp edged entrance of the 
discharge line as a special application of axisymmetric jet 
flow from an orifice, a satisfactory prediction of local 
static pressures along the wall of the discharge line entrance 
was possible. While the jet flow solution was successfully 
used in this one application, the applicability of typical 
jet flow solutions to other geometries could be made which 
would extend the usefulness of the solution presented herein. 

The use of the free streamline solution to obtain a 
modified cavitation number, a*, was shown to have good agree- 
ment with experimental data. The significance of using this 
type of analysis for high velocity flows can be profound in 
application to flows past sharp edged valves of known 
geometries, or past a rounded protuberance extending into 
the flow in a pipe. The free streamline solution predicts 
the values of local velocities which may be very much dif- 
ferent from the typical average value, and which therefore 
require special consideration in any cavitation prediction. 

The development of the time dependent cavitation number 
predicting cavitation in the vibrating discharge line was 
based on a solution of oscillating pressures in a liquid- 
filled flexible tank under no flow conditions. The relations 
for the pressure oscillations were applied to the already 
derived equations for a non-vibrating test section by con- 
sidering the pressure to be a combination of steady plus 
time dependent oscillating components. The velocity was 
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considered to be constant.' This technique of superposition 
was apparent in the final form of the time dependent cavita- 
tion number as it was a combination of the previously obtained 
constant valued modified cavitation number with the addition 
of a time dependent oscillating component. 

The experimental verification of this analytical model 
could not be considered as conclusive evidence to either 
prove or disprove the theory. Since the time dependent cavi- 
tation number was shown to have a sinusoidal form, a phase 
angle relationship between the motion of the test section 
and the strobe lamp which was used to obtain the experimental 
data served as the time reference. Extensive experimental 
work directed towards the time dependency of the cavitation, 
referred to as phase deviation in the text, showed no defin- 
ite trends; however, the predictions for the maximum possible 
variation in the value of the time dependent cavitation 
number were satisfied by the experimental data. 

The lack of phase deviation sensitivity of the experi- 
mental data was attributed to two factors. The first consi- 
deration was the response time of the liquid-vapor evapora- 
tion and condensation processes in the real fluid. The ana- 
lytical model considered a hypothetical liquid which could 
not sustain any tension and had no viscosity; a real liquid 
does not exhibit this type of behavior. The second considera 
tion was the difficulty of visual detection of the instant 
at which cavitation occurred in the liquid. A discernable 
flicker was always noticeable by the observer, and it was 
particularly distracting at low vibrational frequencies. 
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The same free streamline solution used for the cavita- 
tion analysis was used to establish a choking prediction. 

This work , described in Chapter 8, resulted in a determination 
of maximum flow rates which were dependent only on the 
upstream stagnation conditions and the geometry of the 
discharge line. Experimental data correlated well w'ith 
theoretical predictions. 

In order to obtain a relation describing the loss in 
stagnation pressure across the discharge line, it was neces- 
sary to model the separate tw'o-phase flow regimes noted in 
the discharge line. As the liquid entered the sharp edged 
line, the first flow pattern discernable was identical to 
jet flow' from an orifice except that the jet was separated 
from contact with the walls of the line by a vapor annulus. 

At approximately ir-diameters dow r nstream of the entrance, the 
jet flow' underwent a transformation to a two-phase mixture 
most accurately described as being two-phase homogeneous 
flow. The homogeneous regime terminated with a discontinuity 

S' 

where the vapor bubbles suddenly collapsed and formed a liquid. 
The analysis treated this phenomenon as a normal shock front. 
Summing the stagnation pressure losses over the various flow 
regimes yielded an expression for the stagnation pressure loss 
across the discharge line. While the relation was very com- 
prehensive, the void fraction dependency in the discharge 
line did not permit its application to all conditions covered 
in this investigation. However, the special case of minimum 
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pressure loss for, choked flow was shown to agree very well 
with experimental data. 

An expression was developed for the change in stagna- 
tion pressure loss due to a change in the location of the 
normal shock under the assumption of a constant void fraction 
in the homogeneous mixture. When compared to experimental 
data, the solution was shown to be in close agreement. It 
should be noted that the combination of the minimum pressure 
loss relation with the change in pressure loss relation 
yields an analytical prediction for the location of the normal 
shock in the discharge line as a function of the stagnation 
pressure loss on the basis of a constant, previously deter- 
mined, value of the void fraction. While a change in the void 
fraction, was considered to be possible, the magnitude of such 
change was treated as being small enough to permit the assump- 
tion of a constant value. 

While the shock and pressure loss analyses were performed 
only for the situation of flow in the discharge line in the 
absence of vibration, the choking expression for the non-vi- 
brating situation was applied to the experimental data obtain- 
ed from the vibrating discharge line with excellent success. 
Evidentally , the vibration caused no measurable effect on the 
choking of the discharge line. 

The application of the results described in this study 
can be directed towards several areas of technology dealing 
with high speed liquid flows . For example, the cavitation 



parameter formulated in Chapter 7 represents a prediction of 
the limit of applicability of a single phase liquid analysis 
to a high speed liquid flow. The onset of the two-phase 
condition represents a violation of the single phase liquid 
model. As another example, the flashing and subsequent choking 
analyses provided information to predict maximum flow rates 
and minimum pressure drops required to attain these conditions. 
While the relations presented herein are complete and well sub- 
stantiated by an experimental investigation, some specific 
portions of the overall study represent areas where additional 
detailed work can contribute significantly to the understanding 
of the physical processes involved. 

It is recommended that a precise determination of the 
limits of applicability of the jet flow solution used to pre- 
dict the flow pattern at the sharp edged entrance to the dis- 
charge line be made. An experimental study of the pressure 
profiles in the entrance region could be accomplished by means 
of a traversing pressure probe. In this manner, an experi- 
mentally obtained flow pattern could be compared with the flow 
pattern analytically obtained from the free streamline analysis. 

A recommendation for further work on the problem of cavi- 
tation in a vibrating discharge line would be the use of a 
technique for cavitation sensing not directly dependent on 
the physical resources of the observer. A possible facility 
would involve the use of a hydrophone installed very near the 
entrance plane of the line. The signal output could be very 
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carefully filtered so that \vhen the characteristic cavita- 
tion sound was detected, a marker would be signaled to indi 
cate the instantaneous values of pressure and flow rate on 
the' data recorder. This type of detection system would be 
considerably more precise than the visual observation tech- 
nique which was used, especially in determining a phase angle 
relationship . 

Several recommendations can be made to improve the 
understanding of the flow regimes in flashing two-phase floitf 
in the discharge line. The first would involve the use of a 
radiation attenuation technique to experimentally measure 
the void fraction in the discharge line. Photographic tech- 
niques for this purpose are burdensome and inaccurate for 
this type of flow. The void fraction behavior in the immedi- 
ate neighborhood of the normal shock would be most valuable 
as an aid in understanding the physical processes involved. 

The second recommendation deals with the transformation 
of the jet flow pattern at the entrance to the homogeneous 
flow pattern. This could best be accomplished in a two- 
dimensional test section with high speed photographic equip- 
ment. In this manner, perhaps the question of whether nucle- 
ation in the core of the liquid or surface instabilities or 
a combination of both were the principle driving force causing 
the jet deterioration. 

Since the study was concerned only with an adiabatic 
discharge line, the inclusion of heat transfer would 
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significantly add to the usefulness of the analysis. Such 
information as temperature gradients and effects on void 
fraction, pressure losses and flow regimes in the discharge 
line could be investigated not only with water as the working 
fluid but also with some other types of fluids. After a 
thorough understanding of the effects of heat transfer on 
the flow in a non-vibrating discharge line was gained, then 
the extension of the investigation to a vibrating discharge 
line could be made. 
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